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Preface 


This  study  has  been  a very  interesting  and  useful 
learning  experience,  and  I wish  to  convey  my  appreciation 
to  the  Air  Force  Institute  of  Technology  for  requiring  it 
and  to  the  Air  Force  Avionics  Laboratory  for  the  subject 
matter.  I found  that  not  only  is  coherent  combination  of 
lasers  a formidable  engineering  task,  it  is  equally  as  diffi- 
cult to  model  the  problem  on  paper.  This  is  in  part  because 
of  the  cross  talk  between  the  disciplines  of  Electrical 
Engineering  and  Laser  Physics,  but  the  effort  broadened 
(pressure  broadening  I believe)  my  understanding  considerably. 

I had  hoped  to  be  able  to  consider  many  more  applica- 
tions and  in  much  greater  detail,  but  there  was,  as  usual, 
no  time.  I believe  the  results,  however,  to  be  significant, 
particularly  in  the  array  section  where  it  was  found  that  a 
very  large  nvimber  of  lasers  is  needed  to  eliminate  multiple 
lobes  in  the  far  field  combined  beam  pattern.  It  was  also 
found  that,  given  a good  phase  locked  loop,  there  are  few 
major  impediments  to  implementing  a system  of  coherently 
combined  lasers.  This  in  itself  is  a useful  result.  There 
are  however,  many  more  questions  to  be  answered,  and  I hope 
that  this  study  can  be  a useful  starting  point  for  further 
efforts . 

I want  to  express  a great  deal  of  thanks  to  my  advisor. 
Dr.  Stan  Robinson,  for  his  stimulating  discussions.  His 
expertise  and  understanding  of  the  subject  were  invaluable. 
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I also  want  to  thank  my  readers  Dr.  Ted  Luke  and  Dr.  Don 
Shankland  for  their  very  helpful  comments.  Last,  but  not 
least,  I wish  to  thank  my  typists,  Mrs.  Karen  Landis  and 
Mrs.  Rusti  Gaudreau,  who,  though  they  are  still  probably 
wondering  about  all  those  funny  little  symbols,  typed  this 
report  with  great  skill  and  care. 

Hal  E.  Hagemeier 
December  1977 
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Definition  of  Symbols 


A Field  amplitude,  sometimes  indexed,  always 

deterministic 

a,  b Real  and  imaginary  parts  of  q 

B Power  ratio  of  threshold  level  to  theoretical 

average  side  lobe  level 

D,  d Element  spacings  in  an  array,  sometimes  indexed, 

may  be  deterministic  or  random 

E Energy  of  signal  s(t) 

i^y(  ) dimensional  Fourier  transform  on  the  variables 

X,  y 

f^  Frequency  offset  of  ith  laser 

Af  Constant  frequency  offset 

f^  Optical  frequency 

fg  Scan  frequency  (rate) 

g Scale  factor  for  linear  frequency  modulation 

h Peak  spatial  side  lobe  level  with  respect  to  the 

main  beam 

k Propagation  constant,  ^ (k  is  sometimes  used 

as  an  index)  ^ 

L Length  of  laser  array 

m Number  of  pulses 

n Array  parameter 

N Number  of  lasers 

Nq  Twice  the  height  of  the  power  spectral  density 

of  white  noise 

0£  Misalignment  angle  of  optic  axis  for  the  ith  laser 

P(Xj,y,)  Phase  factor  from  the  Huygens-Fresnel  principle 
q Logarithm  (base  10)  of  L in  wavelengths 

q Complex  radius  of  curvature  of  gaussian  spherical 
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R 

R(t) 

s(t) 

S(f) 

tr. 


Uc 

U, 

U, 

u 

w 

w. 


os 


iS 


z 

Z, 


“x»“y 

6 

Y 

6 


e 

C 


Scale  factor  that  determines  the  size  of  the  scan 
sector 

Radius  of  curvature  of  the  phase  front  of  a 
gaussian  spherical  wave 

Correlation  of  a random  process,  may  be  subscripted 
to  indicate  the  particular  process 

Signal  of  interest 

Fourier  transform  of  s(t) 

Temporal  null-to-null  pulse  width 

Temporal  pulse  recurrence  interval 

Coherence  time  of  (P(t) 

i 

Field  at  output  of  laser 
Spatial  part  of  Uq 
Field  at  Zj 
Spatial  part  of  Uj 

Spot  size  of  gaussian  spherical  wave 
Spot  size  at  Zj 

Spatial  null-to-null  pulse  width 

Spatial  period 

Argument  for  z transforms 

Distance  beam  propagates  in  z direction  (no  rela- 
tion to  z transforms) 

Phase  angles 

Confidence  level  of  B 

RIIS  bandwidth  of  s(t) 

RMS  duration  of  s(t)  (no  relation  to  the  Dirac 
delta  function) 

Difference  between  random  variable  and  its  mean 
Beam  angle  in  plane  of  array 
Attenuation  due  to  o 
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e 

X 

V 

T 

4-(t) 

4> 


X 

u 


Phase  angle,  sometimes  indexed,  always  deter- 
ministic 

Wavelength 

Doppler  shift  (frequency) 

Variance  of  a random  quantity  which  is  determined 
by  a subscript,  may  be  indexed 

Time  difference,  t - t^ 

Phase  angle,  always  a random  process,  sometimes 
indexed 

Characteristic  function,  sometimes  subscripted 
Ambiguity  function  (a  two  dimensional  correlation) 
Angle  off  of  array  normal 
Angular  frequency 
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Abstract 


The  coherent  combination  of  several,  single  mode  lasers 
can  produce  a field  similar  to  that  of  a mode- locked  laser 
but  with  more  flexibility.  The  field  for  a quasi-monochro- 
matic  wave  is  considered  a complex,  coherence  separable 
random  process.  The  ensemble  mean  and  covariance  are  deter- 
mined for  the  case  of  a temporally  stabilized  amplitude  and 
a temporal  phase  that  with  the  appropriate  assumptions  is  a 
stationary,  gaussian  random  process.  Mean  fields  are  used 
throughout  as  the  "signals"  of  interest.  The  Huygens-Fresnel 
principle  is  used  to  investigate  the  field  properties  in  the 
Fraunhofer  region  for  two  cases.  The  first  case  is  for  N 
beams  superimposed  with  optical  axes  coincident.  With  appro- 
priate assumptions,  the  performance  of  such  a system  is 
determined  from  the  ambiguity  function.  Range  measurement 
precision  is  found  to  be  proportional  to  where  Af  is 

the  frequency  difference  between  adjacent  lasers.  Velocity 
measurement  precision  is  found  to  be  proportional  to  /NAf^. 
The  second  case  is  for  N beams  in  a linear  array.  The  far 
field  result  is  a scanning  beam  that  in  certain  cases  can 
be  steered.  Aperiodic  arrays  are  considered  in  an  effort  to 
reduce  grating  lobes.  An  array  of  about  1000  lasers  is 
needed  for  reasonably  low  side  lobes. 

Some  generalized  results  are  presented  for  the  effects 
of  misalignment  of  the  lasers.  In  addition,  the  effects  of 
the  random  phases  of  the  lasers  are  considered,  and  the  time- 
space  dual  nature  of  the  results  is  discussed. 
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FIELD  PROPERTIES  OF  MULTIPLE, 

COHERENTLY  COMBINED  IJVSERS 

I.  Introduction 

The  uses  of  laser  systems  in  industry,  medicine,  th^ 
scientific  community,  and  the  military  have  been  growir^g 
steadily  in  the  past  few  years.  Many  applications  have  uti- 
lized lasers  that  are  scanned  or  pointed,  have  high  peak 
powers,  and  are  modulated.  The  military,  for  example,  is 
developing  lasers  for  possible  use  in  ranging  systems,  target 
designators,  terrain  mapping  systems,  and  fire  control  sys- 
tems. One  of  the  types  of  laser  used  belongs  to  the  general 
class  of  mode- locked  lasers. 

Backgroxmd 

Mode- locking  of  an  inhomogeneous ly  broadened  laser  is 
often  done  in  applications  that  require  very  short  pulses 
and  high  peak  powers.  An  inhomogeneous ly  broadened  laser 
typically  oscillates  at  several  longitudinal  modes.  The 
modes  are  at  frequencies  that  are  separated  by  an  amount 
that  is  inversely  proportional  to  the  length  of  the  cavity 
(Ref  33:348-352,365).  Generally,  these  modes  have  random 
phases  with  respect  to  each  other  (relative  phases)  and  the 
resultant  total  laser  output  power  is  proportional  to  the 
number  of  modes.  If,  however,  the  individual  modes  maintain 
a constant  relative  phase,  the  laser  output  power  is 


r 

"i 

[ periodic  with  a period  proportional  to  the  cavity  length,  and 

the  output  peak  power  is  now  proportional  to  the  square  of 
the  number  of  inodes  (Ref  39:256-262).- 

Unfortunately , since  the  pulse  properties  of  the  mode- 
locked  laser  are  determined  essentially  by  the  cavity  length, 
variations  in  pulse  format  are  not  easily  accomplished.  To 
obtain  more  flexibility,  it  seems  reasonable  to  consider 
combining  the  outputs  of  several,  single  mode  lasers.  In 
this  case,  the  phase,  frequency,  and  amplitude  of  each  laser 
can  be  chosen  to  synthesize  desirable  combined  beam  proper- 
ties. If  the  frequency  difference  between  any  two  lasers  is 
an  integer  multiple  of  a constant  frequency  offset  and  if 
the  fields  are  combined  coherently,  the  result  duplicates  the 
output  of  a mode- locked  laser.  Now,  however,  the  frequency 
I offset  can  be  varied  over  a much  larger  range  than  possible 

with  a mode- locked  laser  so  that  the  pulse  properties  of  the 
combined  beam  can  be  chosen  to  meet  more  diverse  require- 
ments.  The  number  of  beams  that  are  combined  can  also  be 
varied  to  change  the  output  pulse  width  and  amplitude.  In 
addition,  the  lasers  could  be  separated  in  space,  i.e.  placed 
in  an  array  of  some  sort,  so  that  spatial  pulse  properties 
exist  analagous  to  the  temporal  properties  of  the  mode- 
locked  laser. 

Analysis  of  the  coherent  combination  problem  has  been 

i 

largely  experimental  and  based  only  on  deterministic  models 

I 

' of  the  individual  laser  outputs  (Ref  1,  16,  17,  41).  There 

j has  also  been  some  speculation  on  the  applicability  of  a 
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"phased  array"  of  lasers  (Ref  7) . Although  these  are  neces- 
sary first  steps,  a more  complete  analytical  approach  which 
considers  the  instabilities  of  the  various  lasers  and  exam- 
ines the  resulting  electromagnetic  field  distribution  in 
the  far  field  is  needed. 

Problem 

This  study  will  focus  on  two  of  the  many  aspects  of 
this  problem.  First,  the  laser  outputs  will  be  modeled  as 
stochastic  processes,  and  the  space- time  electric  field 
distribution  of  the  far  field  will  be  determined.  Second, 
performance  criteria  for  potential  applications  will  be 
determined  and  related  not  only  to  the  characteristics  of 
the  lasers  (such  as  frequency,  phase,  and  spatial  location), 
but  also  to  the  total  nimber  of  lasers.  The  purpose  of  the 
study  is  to  examine  the  effects  of  practical  quantities 
such  as  coherence  time,  phase  instabilities,  number  of 
lasers,  and  position  inaccuracies  on  the  far  field  beam  pat- 
tern and  to  determine  from  a systems  standpoint  what  types 
of  signals  can  then  be  synthesized. 

Scope  and  Assumptions 

Although  most  of  the  calculations  (except  for  specific 
numerical  examples)  will  be  generally  applicable  to  any 
single  mode  laser,  the  underlying  application  will  be  for 
carbon  dioxide  (COj)  lasers.  In  particular,  the  waveguide 
COj  laser  possesses  certain  desirable  characteristics  (Ref 
9) : it  is  homogeneously  broadened  and  therefore  oscillates 


at  one  frequency;  it  is  tvmeable  over  a large  frequency 
bandwidth;  as  a CO^  laser,  it  has  a relatively  long  wave- 
length; and  there  is  a convenient  atmospheric  "window"  at 
the  COj  wavelength  of  10.6  microns.  For  simplicity,  scalar 
fields  will  be  used  instead  of  vector  fields,  and  the  field 
descriptions  will  be  classical  (such  as  the  treatments  found 
in  Ref  12:37-50;  14:11-26;  and  4:10-36).  Furthermore,  the 
Huygens-Fresnel  principle  will  be  used  to  determine  the  far 
field  patterns  (Ref  4:370-375;  12:33-70),  and  a perfect 
channel  (i.e.  free  space)  will,  be  assumed  for  propagation. 

No  attempt  will  be  made  to  consider  the  problems  of  actually 
combining,  phase  locking,  or  controlling  the  phases  of  the 
individual  beams  as  these  topics  are  addressed  in  a related 
s tudy . 

The  desired  result  is  a first  and  second  moment  descrip- 
tion of  the  field  as  a random  process.  The  ensemble  mean 
and  covariance  of  the  field  will  be  considered  as  functions 
of  time  and  space.  The  mean  of  the  field  will  be  used  as 
the  "signal"  of  interest  for  determining  performance  cri- 
teria and  other  results.  The  covariance  will  be  a measure 
of  how  closely  the  actual  results  will  approach  the  mean 
results  and  is  similar  to  the  mutual  coherence  functions 
described  by  other  authors  (Ref  4:499-518;  3:765;  21:235- 
242;  27:359-399;  29:784-789). 

Signal  detection  and  estimation  will  not  be  addressed 
in  detail;  desirable  signal  formats  will  be  those  which  pro- 
duce workable  patterns  within  the  constraints  of  the 
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stochastic  far  field  model.  One  measure  of  performance  will 
be  the  ambiguity  function  (Ref  6:59-105;  31:70-75)  from  which 
the  Cramer-Rao  lower  bounds  on  a given  signal's  measurement 
precision  of  range  and  velocity  can  be  determined.  Another 
perfoimiance  measure  will  be  the  beam  quality  as  a function 
of  laser  position  and  phase  variations.  In  addition,  the 
far  field  beam  pattern  itself  as  a function  of  various  time 
and  space  conditions  will  be  an  indication  of  the  desira- 
bility of  those  conditions  for  certain  applications. 


The  study  is  divided  into  three  sections.  The  first 
section  presents  a statistical  model  of  the  individual  laser 
modes.  They  are  considered  to  be  complex,  coherence  separ- 
able processes  whose  real  and  imaginary  parts  are  identically 
distributed.  For  practical  cases,  the  mode  amplitudes  are 
considered  deterministic,  and  the  randomness  (laser  insta- 
bilities) appears  in  the  phase  of  each  beam.  The  duality  of 
time  and  space  variations  is  shown  as  a general  case  when 
several  beams  are  combined.  Far  field  results  are  deter- 
mined and  the  effects  of  physical  limitations  are  studied. 
Finally,  two  general  classes  of  combined  beams  are  consid- 
ered: the  coaxial  case,  and  the  side-by-side  array  case. 

The  means  and  covariances  are  determined  for  both  cases, 
and  some  characteristics  of  the  means  are  examined. 

The  second  section  considers  the  coaxial  case  in  more 
detail.  Coherent  detection  and  a matched  filter  receiver 
are  assumed,  and  the  ambiguity  function  of  the  combined 
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field  mean  is  determined.  This  shows  explicitly  the  depend- 
ence of  the  measurement  precision  of  the  signal  as  a function 
of  the  nvimber  of  lasers,  the  frequency  differences  between 
lasers,  and  the  relative  phases  of  the  lasers.  Potential 
applications  for  ranging,  velocity  measurement,  and  both  are 
considered  in  light  of  the  performance  criteria. 

The  third  section  presents  analytical  results  for  var- 
ious linear  arrays.  The  scanning  characteristics  of  the 
combined  beam  in  the  far  field  are  studied.  Pattern  synthe- 
sis and  other  beam  properties  are  related  to  the  number  of 
lasers  and  their  position.  In  addition,  aperiodic  arrays 
and  the  statistics  of  random  arrays  are  included  in  an 
attempt  to  reduce  high  side  lobes. 
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II.  Field  Models 

The  optical  fields  used  here  will  be  described 
classically.  Most  generally,  the  field  for  a monochromatic 
linearly  polarized  wave  may  be  written  as 

u(x,y,z,t)  = A(x,y,z,t)cos(2TTf5t-4)(x,y,z,t)) 

where  u(x,y,z,t)  is  a scalar  function  representing  either  an 
electric  or  a magnetic  field  at  position  x,  y,  and  z,  and 
time  t;  A is  the  amplitude  of  the  wave;  <J)  is  the  phase;  and 
f,  is  the  optical  frequency.  Typically,  u is  written 

u(x,y,z,t)  = Re[U(x,y,z,t)exp(-j2Trf5t)) 

where  Re  means  "the  real  part  of"  and  U(x,y,z,t)  is  a complex 
function  called  the  complex  envelope  of  u(x,y,z,t),  V7here 

U(x,y,z,t)  = A(x,y,z,t)exp(j(})(x,y,z,t)) 

It  is  assumed  that  U is  a complex  random  process  whose  real 
and  imaginary  parts  are  identically  distributed  and  that  z 
is  the  direction  of  propagation.  The  field,  then,  at  any  z 
is  a function  of  x,  y,  and  t and  is  designated  by  a numerical 
subscript  on  U,  for  example 

U(x,  y,  Zj,t)  = U,(x,y,t) 

Output  Field  Model 

U can  now  be  written 

U(x,y,t)  = A(x,y,t)exp(j(t)(x,y,t)) 
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v/here  the  amplitude  and  phase  are  slowly  varying  random 
functions  of  time  and  space.  For  a single  mode  laser,  the 
amplitude  is  usually  considered  spatially  coherent  for  a 
given  z and  separable  (Ref  3:768,770)  so  that 

A(x,y,t)  = Ag(x,y)A(t) 

where  Ag  is  the  spatial  dependence  of  the  field  amplitude 
and  is  generally  considered  to  be  knovm.  (The  last  section 
will  deal  with  results  for  x and  y unknown.)  A(t)  is  the 
randomness  of  the  amplitude  due  mostly  to  spontaneous  emis- 
sion (Ref  3:765;  21:236).  If  the  laser  output  is  considered 
to  be  amplitude  stabilized  (Ref  2:32;  18:378;  28:2399;  29:784; 
and  33:63),  then  A(t)  becomes  a simple,  non-random  constant, 
A,  and  the  fluctuations  of  the  laser  light  are  primarily  in 
phase  (the  effect  of  a random  Ag(x,y)  will  be  considered 
later).  In  addition,  (p(x,y,t)  is  also  considered  separable 
so  that 

<()(x,y,t)  = (t>g(x,y)  + (j)(t) 

where  ((>  is  the  non-random  spatial  dependence  of  the  phase  of 
s 

the  field  (the  random  case  is  treated  later)  and  <|)(t)  is  the 
randomness  of  the  phase  that  accounts  for  the  line  broadening 
of  the  laser.  Therefore, 

U(x,y,t)  = A{As(x,y)exp(j(^g(x,y))  }exp(j4.(t)) 

= A U (x,y)  exp(j(f(t))  (1) 

s 
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It  is  assumed  that  the  single  mode  lasers  discussed  here 
operate  in  the  TEM^^  mode  so  that  the  spatial  part  of  the 
field,  Ug(x,y),  has  the  form  of  a gaussian  spherical  wave, 


UgCx.y) 


(2) 
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where — is  a normalizing  factor  for  unit  power  flow  (i.e. 
00 

// lUg I ^dxdy=l) , R is  the  radius  of  the  spherical  wave  phase 

front,  and  w is  the  "spot  size"  of  the  gaussian  amplitude 
distribution  (equivalent  to  the  e“^  point  or  /2  times  the 
standard  deviation  of  a gaussian  curve)  (Ref  33:307). 

Single  Beam.  From  Eqs  (1)  and  (2)  for  the  output  of  a 
single  laser  at  some  initial  plane  in  front  of  the  laser 
output  mirror, 


Uo(x,y,t)  Aexp(j  (e+(|)(t))  (3) 


where  q is  the  complex  radius  of  curvature  defined  by 


1=1- 
q R '^irw 


(4) 


and  0 has  been  added  as  an  arbitrary  phase  angle.  This  phase 
angle  is  considered  a controllable  parameter  (it  could  be  a 
function  of  time)  which  will  prove  useful  later  on.  A is 
also  assumed  to  be  a controllable  quantity. 

The  random  phase  of  a free-running  laser  is  ultimately 


I 


I 


9 
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i 

a result  of  spontaneous  emission  and  quantvim  noise.  These 
noise  sources  cause  the  phase  to  be  a stochastic  process 
known  as  a random  walk  (Ref  18:370-371;  34:181).  It  is  well 
known  that  such  a process  can  be  described  as  a Wiener-Levy 
process  (sometimes  called  Brownian  motion)  that  is  zero  mean, 
gauss ian,  and  non-stationary . Therefore,  using  known  results 
for  the  characteristic  function  of  a gauss ian  random  variable 
(Ref  26:254,474-510), 

E(Uj(x,y,t))  = AUjs(x,y)E  exp{  j (0+({)  (t))  } 

= AU5g(x,y)exp  je  - — (5) 

where  E(  ) means  to  compute  the  ensemble  mean  or  expected 
value  of  the  quantity  in  brackets,  and  a*  (t)  is  the  variance 
of  the  random  process  (()(t).  If  n(t)  is  a white  noise  pro- 
cess that  is  zero  mean,  gaussian,  and  stationary  then  the 
Weiner  process,  (J)(t)  can  be  written  (Ref  26:502), 

<l>(t)  = n(u)du 

If  a is  the  height  of  the  power  spectral  density  of  n(t)  , 
then  the  variance  of  (j)(t)  is, 

a*(t)  = at  (5a) 

Therefore,  the  mean  of  in  Eq  (5)  goes  to  zero  as  t 
increases. 
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The  correlation  of  now  becomes 

E(Uo(x.y.t)U*(x^y^t'•)) 

» * 

=E  AU5g(x,y)exp{j  (0+<j)(t))  }AU*g(x'',y')exp{-j  (e+(J)(tO)  } 

=A2U5g(x,y)U*g(x'',y')E  exp{j  (({) (t) -cj) (t") ) } 

V • 

=A^U,g(x.y)U*g(x^yOexp[- 

=A2Uos(x,y)U^s(x",y")exp  - ^|t-t"|  (6) 

where  o\  is  the  variance  of  the  gauss ian  random  process 

A(},(t.tO  = <}>(t)  - (7) 

-2 

and  - ^ (Ref  8:139;  26:476,  502).  Now  the 

correlation  of  U,  in  Eq  (6)  is  stationary  even  though  the 
phase  (J)(t)  itself  is  not. 

Multi-Beam  (Coaxial) . Figure  1 (on  the  next  page)  shows 
the  case  where  N beams  are  combined  coaxially.  No  consider- 
ation is  given  here  to  how  the  beams  are  combined.  Since 
this  investigation  was  motivated  by  a mode-locked  laser,  each 
mode  can  be  modeled  by  a single  mode  laser  with  a frequency 
offset,  so  that  a frequency  term  is  added  to  the  time  depend- 
ent part  of  Eq  (3) , 

Uj(x.y,t)  = U5g(x,y)  Ai  exp{ j (2Trfit+e^+4>^(t))  ) (8) 
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Individual 
Output  Beams 


Combined  Output 
Beam 


(x,y,t) 


Fig.  1.  Multi-Beams  Combined  Coaxially 

where  is  the  (controllable)  frequency  spacing  between 

adjacent  spectral  components  of  Uj.  The  number  of  lasers,  N, 

N-1 

that  are  used  is  only  limited  by  the  total  bandwidth,  f^, 

since  the  total  bandwidth  must  be  less  than  or  equal  to  the 
gain  line  width  of  the  lasers  being  used.  In  general,  f^ 
could  have  an  additive  random  comp.onent,  but  this  will  be 
ignored  here  for  simplicity.  It  is  also  assumed  that  the 
♦ i(t)  are  statistically  independent.  Therefore,  as  with 
Eqs  (3)  and  (5) , 


E(U,(x,y,t))  = U.g(x,y)  Z Ai  exp |j  (2TTf ^t+Bi)  - (9) 


'I 


where  a?(t)=a^t  is  the  variance  of  the  phase  of  the  ith  beam. 
This  might  be  called  an  "open  loop"  expression  (similar  to 
Eq  (5))  in  that  Aj^,  fj_,  and  0j^  are  deterministic  and  unre- 
lated to  i.e.  there  is  no  attempt  to  use  measurement 

dependent  feedback  to  compensate  for  <j)jj^(t).  The  covaria  ce 
of  Uj  is 

Cov(Uo.Uj} 

= E(U„(x,y.t)U*(x^y^tO)  - E (U,  (x.y . t)]  E (U*(x^y  ^ t ') ) 


= U5s(x,y)U*g(x",y")  Z A?  exp  (j 2Trf i(t- 1") ) 


{exp 


W • I 


exp 


a. 


(10) 


Although  the  covariance  is  non- stationary  (since  the  mean  of 
(|)(t)  is  non-stationary)  , the  last  term  in  Eq  (10)  goes  to 
zero  quickly  with  t and  t'  and  the  covariance  is  quasi- 
stationary. The  mean  of  the  field  in  Eq  (9)  behaves  as  in  Eq 
(5)  , and  the  covariance  of  in  Eq  (10)  varies  with  the 
statistics  of  <|)^(t)  similar  to  Eq  (6). 

It  is  possible  to  continue  using  the  non-stationary 
statistics  of  the  free-running  phase  (l)j|^(t)  as  is  done  in 
Eqs  (9)  and  (10).  However,  this  will  render  the  results 
unusable  for  certain  applications  (notably  the  ambiguity 
function  in  the  next  section).  In  addition,  there  is  pres- 
ently no  known  way  to  measure  such  field  quantities  directly. 
Therefore,  the  following  assumptions  will  be  made:  First, 
the  mean  of  the  field  (as  opposed  to  the  field  intensity) 
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will  be  considered  the  signal  of  interest.  Second,  it  will 
be  required  that  for  the  field  of  any  laser  to  be  observed 
at  any  point  in  space  or  time,  it  must  be  coherently  detected. 
This  necessitates  a receiver  structure  as  shown  in  Fig.  2. 

In  general,  with  no  control,  the  phase  of  the  detected  signal 
would  be  proportional  to  (J>^(t)  - (})j(t).  Since  both  are 
independent  Wiener  processes,  their  difference  is  also  a 
Wiener  process.  However,  if  the  control  is  used,  then  the 
local  oscillator  phase  (ti,(t)  becomes  a function  of  a measure- 
ment made  by  the  closed  loop  receiver.  This  measurement  is, 
in  turn,  a function  of  the  phase  difference  <J)^(t)  - (j)j(t). 
Therefore,  <|>g(t)  is  no  longer  independent  of  (J)j^(t).  If  the 


Fig.  2.  Coherent  (Heterodyne)  Detection  of  a Field 


i 

i 

j control  apparatus  is  properly  designed,  the  phase  difference 

^ 3 stationary  process  since  the  control 

mechanism  is  tracking  the  temporal  variations  of  the  differ- 
ential phases,  i.e.  4>p(t)  is  nearly  locked  to  <^^(t).  (This 
assumption  is  often  made  implicity;  see,  for  example.  Ref  21; 
237;  24:78,85.)  In  an  absolute  sense  the  two  phases  are 
still  proceeding  in  a random  vjalk  manner,  however,  the  receiv- 
er is  only  observing  the  difference  betvzeen  the  two  phases, 
and  this  difference  can  be  controlled  arbitrarily  v/ell  in 
theory.  It  is  not  generally  possible  to  make  precisely 

equal  to  <j)^(t)  thereby  making  the  observed  signal  totally 
non-random.  Thus,  the  receiver  of  Fig.  2 usually  produces 
non-exact  or  noisy  phase  control,  wherein  some,  but  not  all, 
of  the  received  field’s  instabilities  are  compensated. 

For  the  multi-beam  case  in  Eq  (8) , when  the  field  is  to 
be  observed,  a receiver  such  as  the  one  illustrated  by  Fig.  2 
must  be  used  to  coherently  detect  the  signal.  To  simplify 
matters  Eq  (8)  can  be  rewritten  in  terms  of  measureable  field 
quantities  as 

Uj(x,y,t) 

= Ujg  (x.y)  Ai  exp{j  (2TTfit+e^+4»^(t)-(|)5  (t))  } 

= Ujg  (x,y)  Ai  exp{j  (2ufit+e^+(t't (t)) } (11) 
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where  <!)r(t)  is  defined  to  be  the  difference  between  the  phase 
of  the  reference  laser,  ())o(t),  and  the  phase  of  the  ith  laser. 
If  the  phase  of  each  laser  is  locked  to  the  reference  (simi- 
lar to  Fig.  2)  then  each  of  the  in  Eq  (11)  can  be 

considered  stationary.  The  use  of  a phase  locked  loop  as  a 
control  mechanism  will  invariably  change  the  gauss ian  random 
process  (})£(t)  into  a non-gaussian  process  (J'^(t)  such  that 

E (exp{  j (|)C  (t)  })  is  not  necessarily  exp  - — ^ and  the  (|)^(t) 

are  not  necessarily  independent.  Nevertheless,  for  the  case 
of  a strong  signal,  i.e.  for  a phase  locked  loop  with  a high 
signal- to-noise  ratio,  the  (p^(t)  may  be  modeled  as  independ- 
ent gaussian  random  processes  (Ref  38:86-92). 

Since  all  of  the  lasers  are  now  locked  temporally  to  the 
phase  of  the  local  oscillator,  the  phase  of  the  signal  of 
interest  becomes  the  controlled  phase  differences  between  the 
individual  lasers  and  the  local  oscillator,  and  it  is  sta- 
tionary. Hence  the  random  walk  nature  of  the  phase  of  a 
free-running  laser  (or  lasers)  will  always  be  transparent  as 


far  as  any  results  derived  from  Eqs  (8)  or  (11)  are  concerned. 
By  definition,  then,  Eq  (8)  will  be  understood  to  have  4)^  (t) 
(=(|)^(t))  that  are  zero  mean,  stationary,  independent,  gaussian 
random  processes  when  observed  in  the  required  manner  (Fig. 2). 
The  mean  of  the  field  now  becomes 

E(Uo(x,y,t))  = Vg^(x,y)  I Ai  exp |j  {2TTfit-bO^}  - (12) 
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where  the  variance,  a?,  is  a constant,  and  the  covariance 
of  U.  becomes 

Cov(U5,Uj)  = Ujg(x,y)U*g(x",y'')  I A?  exp  (j2TTfi(t-t")) 

{exp(-a?)  (exp(R({>^(T))-l)}  (13) 

where  R((ij^(T)  is  the  correlation  function  of  ij)j^(t)  and  depends 
only  on  t = t-t'.  Now,  o?  is  a measure  of  how  well  the  ith 
laser's  phase  is  controlled.  The  degree  of  control  necessary 
can  be  calculated  based  on  system  requirements  (this  will  be 
done  later) . It  should  be  noted  that  although  Eqs  (12)  and 
(13)  are  apparently  statistics  of  the  field  of  interest,  they 
are  technically  -the  statistics  of  the  product  of  the  field  of 
interest  and  the  field  of  the  local  oscillator  in  the  required 
receiver.  Since  this  assumption  is  always  made,  however, 
explicit  use  of  the  local  oscillator  field  will  not  be  used. 

Multi-Beam  (Array).  It  may  not  be  possible  or  desirable 
for  physical  or  economic  reasons  to  combine  the  individual 
laser  beams  coaxially.  In  this  case  the  lasers  may  be 
aligned  side-by-side  in  a one  or  two  dimensional  array  simi- 
lar to  the  antennas  in  a phased  array  radar  (Ref  7:1103). 
Figure  3 (on  the  next  page)  illustrates  a rectangular  array 
where  is  the  (center- to-center)  spacing  in  the  x direc- 
tion, d^j^  is  the  spacing  in  the  y direction,  and  and  d^^ 
are  larger  than  the  spot  size,  w,  for  all  i and  k.  Now,  as 
. shown  in  Eq  (14)  on  the  next  page,  the  field  becomes 

I 

i 
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Fig.  3.  Front  View  of  Laser  Array 


n-1  m-1 


' ifo  k=0  = ''ik 


exp{j  ) (14) 

where  mn  = N.  The  lasers  could  be  arranged  in  a circle, 
spiral,  or  any  other  geometrical  shape,  but  for  simplicity 
the  array  is  considered  rectangular  unless  otherwise  noted. 
Essentially  Eq  (14)  is  the  same  as  Eq  (8)  except  that  the 
spatial  term  is  included  within  the  summation.  Similarly, 
the  mean  of  the  field  is. 
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and  the  covariance  of  Uj  is 


Cov(U.,u;)  - E E ^ik 

{exp(j2itfj^^(t-t')){exp(-oU  (exp(R4.j^j^(T))-l))  (16) 

Since  the  spatial  and  temporal  parts  of  the  field  have  been 
assumed  to  be  separate,  the  mean  and  covariance  of  the  field 
in  Eqs  (15)  and  (16)  behave  the  same  with  respect  to 
as  Eqs  (12)  and  (13),  i.e.  only  the  spatial  part  of  the  field 
has  been  changed. 

Far  Field  Model 

The  fields  given  in  Eqs  (3) , (8) , and  (14)  are  now  prop- 
agated through  space  by  means  of  the  well-known  Huygens - 
Fresnel  principle  (Ref  4:370-375;  12:40-46).  Interactions  of 
the  field  with  the  atmosphere  are  not  considered.  In  addi- 
tion, the  Fraunhofer  condition  (Ref  4:382-386;  12:61)  is 
imposed  since  consideration  is  given  here  only  for  the  far 
• field.  (Since  the  space  and  time  parts  of  the  fields  are 

separate  and  the  Huygens-Fresnel  principle  affects  only  the 
spatial  function,  the  temporal  results  that  will  be  deter- 
j mined,  including  the  effects  of  the  random  process  (f^(t), 

I are  unchanged  by  propagation  through  space.  If  desired,  the 

I 
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fields  could  be  determined  in  the  near  field  using  the 
Fresnel  condition  as  in  Ref  4:382-383.  However,  v/ith  little 
loss  of  generality,  the  Fraunhofer  condition  is  used  here  for 
simplicity.)  A sufficient  condition  for  the  Fraunhofer  case 
requires  that 


z,»- 


irCx^+y^) 


This  may  be  a rather  stringent  condition,  for  example,  in 
the  multi-beam  array  case  where  (x+D.,  and  (y+d., 

may  be  on  the  order  of  10  cm.  Then  if  X = 10. by,  must 
be  about  600  m to  be  an  order  of  magnitude  larger  than 


After  propagating  a distance  Zj  down  the  optical  axis, 
the  field  becomes 


Uj(Xj,y,,t)  = exp  jkZj  exp  {7  U5(x,y,t) 


2t\ 

exp  -j  (r— )(xx  +yy  ) dxdy 

. AZ  J 


= PCXj.yj)  ^„(U5(x,y,t)l 


where  P(x^,y^)  is  just  the  phase  factor-in  front  of  the  inte- 
gral and  ) is  the  two  dimensional  (spatial)  Fourier 

xy 

transform  of  Uj(x,y,t)  on  x and  y (Ref  12:61).  Although  Uj 
is  considered  a function  of  the  rectangular  coordinates 
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Xj  and  yj,  it  is  also  a function  of  the  spatial  frequencies 

“ (Xi)/(XZj)  and  f^  = (yi)/(AZj)  in 
keeping  with  the  nature  of  the  Fourier  transform.  In  the 
following  discussions,  Uj(x,y,t)  is  not  considered  to  have 
an  explicit  aperture  function  limiting  its  spatial  extent  at 
the  laser  output.  As  noted  in  Ref  33:312-313,  a conservative 
aperture  design  for  maximum  power  transmission  makes  the 
aperture  size  about  3w  (where  w is  the  spot  size  of  the  beam) 

The  diffraction  effects  of  this  aperture  on  the  beam  profile 
are  completely  negligible,  hence  there  is  no  explicit  depend- 
ence of  Uj  on  the  aperture.  Any  other  desired  aperture 
function  may  be  used  at\d  would  become  simply  (though  per- 
haps tediously)  the  spatial  convolution  of  the  Fourier  trans- 
forms of  Ujg(x,y)  and  the  desired  aperture. 

Single  Beam.  If  the  field  of  Eq  (3)  is  used  in  Eq  (18) , 
the  resulting  far  field  for  a single  beam  is 

Ui(x,,yi,t)  = P(Xj,yj)  ,^y(Upg(x,y))A 

exp{j  (0+(f)(t)) } (19) 

The  mean  of  the  far  field  becomes 

1 

I 

E(U,  (Xj  ,yi  ,t))  = P(xi  ,yj)  (Ujg(x,y))  A exp  je-^  (20)  I 

and  the  correlation  becomes 


E(U,(Xj,yj,t)U*(x;,y;,t"))  = P(x, ,yj)P*(x^ ,y;)Uj3(Xj ,y,) 

U*g(xJ,yO  A*  exp  (-o^+R(|)(t))  (21) 


where  UigCxj.yi)  = (Up  g (x,y))  , For  completeness 


> TT 


2 1 


- ^ 
r w '^x 


exp(-j-^  x^) 


\ ^ r 


exp(-j—  y*) 

Xq 


■4^  lif 


(22) 


From  Eq  (20) , 


P(Xj,y,)^j.y(Uos(x.y)) 


--/I 


2 q exp(jkzj 


. > ir  V7Z , 


exp 


(23) 


and  from  Eq  (21) , 


P(Xj  ,yj)P*(x;‘ ,ypUjs(Xj  ,y j )U*s (x^ ,y ^ 


“ rw^  exp{j  Y^((q+Zj)  (x2+yp-(q*+Zj)  (x^2+y-2)j  j (24) 


Xz 


The  intensity  of  this  single  beam  in  the  far  field  is 


I,  = Uj (Xj .yj ,t)U*(Xj ,yj ,t) 


exp 


TTW^zf 


2TTb  /.,2  i„2 


I XZ^ 


2 (x^+yp 


(25) 


where  from  Eq  (4)  ^ has  the  form  a + jb  where 


J 
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r 


( 

k 


Zj  (trw^ ) ^ ^ _ ttACwZ))^ 

(tiw2)2+(az^)2  ' “ (Trw2)2  + (XZj)2 

I and  in  the  far  field  case  R = z^.  In  Eq  (25)  the  far  field 

] intensity  of  the  single  beam  is  completely  deterministic 

since  so  far  the  randomness  has  been  assumed  to  be  only  in 
phase. 

Multi-Beam  (Coaxial) . Again,  the  extension  can  be  made 
to  svim  N beams  and  the  far  field  now  shows  some  interesting 
properties.  Uj(Xj,yj,t)  in  this  case  can  be  determined 
either  by  adding  N versions  of  Eq  (19)  with  the  appropriate 
I addition  of  a frequency  offset,  f^j^,  or  by  applying  the 

i 

j Huygens-Fresnel  principle  to  Eq  (8).  The  same  result  is 

I 

j obtained  in  either  case, 

i 

I 

Ui(Xj,yj,t)  = P(Xj ,y,)Ujg(Xj ,yj)  Z A. 

i 

exp{j(2TTf^t+0j.+(^^(t)]}  (26) 

The  mean  of  this  far  field  is 

E(Uj(Xj ,yj ,t))  = P(Xj ,yj)Uig(Xj ,yj)  Z A^ 

exp[j(2iTf^t+0^)-  ^ (27) 

The  covariance  is 

Cov(U,,ur)  - P(x,.y,)P*(xr,y;)U,5(x.,y,)U*5(x;.y;)  r A| 

exp  {exp (-op  (exp  (R(>^(t))  -1)1  (28) 
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The  intensity  of  this  far  field  is  not  deterministic  as  it 
was  for  a single  beam  in  Eq  (25)  because  of  the  cross  terms 
between  beams.  The  mean  of  the  intensity  is 


E(li)  = 


2(a'+b2) 
-"■-2-1"  exp 
TTW  Z, 


2Trb 

Xz? 


(xf+y^) 


E + ^E^  ^E^  A.  A, 

i 3.  i=0  k=0  1 ^ 

i+k 


exp{j(27r(f.-fj^)t+(e.-ej^))  - ^■■±■^-1^} 


(29) 


This  is  a somev?hat  formidable  expression,  but  some  minor 
simplifying  assumptions  later  will  make  it  more  useful. 

Examination  of  the  mean  in  Eq  (27)  yields  some  interest- 
ing results.  It  is  reasonable  to  assume  that  although  the 
(()^(t)  are  independent,  the  a?  are  certainly  similar  if  not 
identical.  This  is  because  of  the  nature  of  the  noise 
sources  producing  the  random  phase  fluctuations,  because  of 
the  similarity  of  the  phase  locked  loops,  and  because  all  of 
the  lasers  are  assumed  to  be  of  the  same  type  and  construc- 
tion. Also,  since  the  A^  are  controllable,  they  are  assumed 
to  be  the  same  for  convenience.  Therefore,  Eq  (27)  becomes 

E(Uj (Xj ,yj ,t))  = P(Xj ,yj)Ujg(Xj ,yj)A  exp 

E exp  (j  (2TTf . t+e  . ))  (30) 

i ^ ^ 

At  this  point  some  assumptions  must  be  made  for  f^  and 
0^  so  that  a closed  form  solution  can  be  obtained.  If  the 
frequencies  between  the  lasers  are  all  the  same  (as  in  a 


2A 


mode-locked  laser)  then  = iAf  where  A£  is  the  constant 
frequency  spacing.  A similar  assumption  can  be  made  for  0^^ 
but  the  results  will  not  be  affected  in  any  significant  or 
useful  manner,  so  0^  is  considered  zero.  Now,  Eq  (30) 
becomes 

E[Uj (Xj ,yj ,t))  = P(Xj,yi)U  (x,,yj)  A 


« o*  , . rN-lv.,  sin  NirAft  /o- 

exp  - — + j(— )2TTAft  — ^ (3. 

I ^ ^ J sin  irAft 


The  (sin  Nx)/(sin  x)  term  in  Eq  (31)  (the  typical  result  of 
a mode-locked  laser  as  in  Ref  39:131)  is  a pulse  train  in 
time  as  shown  in  Fig.  4.  The  pulses  have  a null-to-null 


sin(a)N/2) 


Fig.  4 Plot  of  sin(N  2) /sin (2)  for  N=8,  u)=2TrAft 
pulse  width  of 

t » — ^ 
p NAf 

and  a pulse  recurrence  interval  (pulse  rate)  of 


For  the  waveguide  CO2  laser  (whose  gain  line  width  is  about 
700  MHz),  NAf  must  be  less  than  or  equal  to  700  MHZ.  There- 
fore the  minimum  achievable  pulse  width,  t^,  is  about  3 nsec. 
Tp  depends  on  Af  which  is  limited  by  the  total  niamber  of 
lasers  used.  The  peak  amplitude  of  the  (sin  Nx)/(sin  x) 
function  at  each  T^  is  just  N as  expected  for  coherent  super- 
position of  N sources.  The  pulse  V7idth  and  pulse  rate  both 
depend  on  N and  Af  which  are  assumed  to  be  under  the  opera- 
tor's control.  Eq  (31)  also  shoves  that  the  apparent  frequency 
of  the  field  is  changed  to  the  optical  frequency  f^  minus 

(^^~)Af.  Ref  41  reports  experimental  results  of  Eq  (31). 

If  the  assumptions  for  Eq  (31)  are  applied  to  the  inten- 
sity then  the  series  terms  in  Eq  (29)  become 


E A?  = A^N 
i 


W exp{j  (2TT(f^-fj^)t+(0^-e^))- 

= A*  exp(-a^)  E E exp  (j  (i-k)  2TrAf t) 
i=|=k 

= A*  exp  (-o^)  {E  (N-i)  (exp  (j  2TrAf  ti)+exp  (- j 2TTAfti)  ) } 


= A^  exp(-o*){2  E (N-i)cos  (i(2TTAft))  } 
i 

.2  f 2'\fl-cos  N2TrAft  „ 

= expl-o-'J  — 5 N 

^ I 2 sin^  TTAft 

- A'  o,xp(-„  = )((2i2-HlM!l)'  nI 

‘ ^ M sin  TTAft 
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(See  Ref  22:104-105  for  series  sums.)  Therefore, 


E(I,)  - 


2(aHb^) 

TTW^Zj 


exp 


2TTb 

'Az^ 


(xj+y,^) 


(l-exp  (-»=))«+  exp 


(34) 


If  exp(-a^)->-0  and  E(Ij)  becomes  a constant  with  respect 

to  time.  This  is  equivalent  to  the  result  of  a multi-mode 
laser  with  no  mode- locking.  If  a^->0,  exp[-0^)^l  and  E(lj) 
is  a pulse  train,  i.e.  E [l  j) ->  1_E  [U^  (x^  ,y  j,t)]  [ ^ for  a^->-0. 

This  is  equivalent  to  a mode- locked  laser  output. 

The  random  phase  fluctuations  <}>^(t)  simply  cause  the 
amplitude  of  the  mean  of  the  field  in  Eq  (31)  to  be  reduced 
by  exp(-  ^ )•  As  the  variations  in  (J)^(t)  become  larger,  the 

amplitude  of  the  pulses  in  the  far  field  becomes  smaller. 

From  this  relation,  the  degree  of  control  of  the  phases  of 
the  lasers  can  be  calculated.  If  the  desired  maximum  atten- 
uation is  n,  then 


^ “ *^"2  In  n 


(35) 


For  example,  if  n is  desired  to  be  no  larger  than  1 db  (about 
0.80)  then  a must  be  about  0.67  rad.  Similarly,  for  phase 
control  on  the  order  of  1 rad,  q is  about  50%.  The  , effect  of 
([)^(t)  becomes  even  clearer  upon  examination  of  the  total 
power  flow  through  the  z = Zj  plane.  The  mean  power  flow  is 


E(Power)  = E // | Uj  (x,  , y j , t)  | Mxj  dy  j ] 

V-00  J 
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E (Power)  = //  E(lj)dXjdyj 
— 00 

Since  Ujg(x,y)  is  already  normalized  for  unit  pov;er  flow  (Eq 
(2))  and  the  Huygens-Fresnel  principle  merely  changes  the 
distribution  of  the  power,  not  the  amount,  the  mean  power  is 

E (Power)  = | (1-exp  (-0^)  )N+exp  (-0^) 

Therefore,  at  t = 0 (pulse  peak)  and  for  A = 1, 

E(Power)  = (1-exp  (-0^)) N + exp(-0^)N^  (36) 

As  o^-*-oo,  E (^wer)  approaches  N (incoherent  superposition)  and 
as  0*-»-O,  E (Power)  approaches  (coherent  superposition). 

It  is  also  instructive  to  examine  the  behavior  of  the 
covariance  of  Uj  in  Eq  (28),  For  0^^  = 0, 

Cov(Uj,Ui)  « exp  (-0^)  (exp  (R(j)(T)) -1)  (37) 

The  general  characteristics  of  R<|)(t)  are  that  R(f)(0)  = o^  and 
that  R4>(t)  = R(|)(-t),  R(() (0) ^ I R<j) (t)  I , and  generally  R())(t)-^0 
as  T-»^  (since  <(>jj^(t)  are  zero  mean  processes)  (Ref  23:126; 
26:336-338).  At  t = 0,  the  covariance  of  becomes  the 
variance,  c^,  of  and 

o,*  « 1-exp  (-0^)  (38) 

By  the  Cht5byshev  inequality  (Ref  26:150)  the  probability  that 
a random  variable  differs  from  its  mean  by  more  than  an 
amount  e is  less  than  or  equal  to  its  variance  divided  by  e*. 
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Therefore, 


Pr(|U,-E(uJ  |>c)<  ^ 

Now  as  o^-»-0,  Pr  ( I Uj -E  (Uj ) I ^e)  approaches  zero,  and  Uj  is  equal 
to  its  mean  with  probability  one,  i.e.  the  smaller  the  vari- 
ance of  <ti^(t)>the  more  closely  the  field  approaches  its  mean 
in  Eq  (27) . 

As  T->“,  the  two  sample  functions  (f)(t)  and  (f)(t')  become 
more  and  more  uncorrelated.  Typically,  R(i)(T)  = 0 for  all 
t^T^  where  T^  is  called  the  coherence  time  of  the  process 
(|)(t).  For  time  differences  x larger  than  the  coherence 
time,  the  phases  at  the  two  times  are  uncorrelated.  From  Eq 
(37) , t>T^  implies  that 

Cov(Uj,Ui]  = 0 

i.e.  the  coherence  time  of  the  phase  <|>(t)  is  the  coherence 
time  of  the  field  Uj.  In  general,  Cov[Uj,Uj)  will  decrease 
faster  with  x than  will  R(})(x)  so  that  the  effective  coherence 
time  of  Uj  may  be  less  than  that  for  4)(t),  however,  this 
is  less  true  as  o^-*Q  (Ref  13:66).  The  coherence  length  of 
the  field  is  related  to  the  coherence  time  by  the  speed  of 
light.  (A  longer  coherence  time  implies  a greater  coherence 
distance.)  Therefore,  field  points  separated  by  tim.es  (dis- 
tances) greater  than  the  coherence  time  (distance)  are 
uncorrelated. 

One  further  comment  on  Eq  (27)  concerns  the  beam  waist. 
Since  the  individual  beams  are  added  coaxially,  the  waist  of 


the  total  beam  propagates  as  the  waist  of  any  single  beam. 
Therefore,  the  waist  Wj  at  a distance  is 


Jill 

ttw 


(39) 


where  w is  the  waist  at  z = z,  (Eq  (2)),  and  the  far  field 

2 

here  implies  (Ref  33:308).  For  example,  if  w = 1mm 

A 

and  X = 10. 6u,  then  z>>0.3  m which  is  a much  less  stringent 
condition  than  the  Fraunhofer  condition  of  Eq  (17) . There- 
fore, Eq  (39)  is  essentially  an  equality  when  the  Fraunhofer 
condition  is  met. 

Multi-Beam  (Array) . For  this  case,  Eq  (26)  becomes 


U,(x,,y,,t)  = I Z P(x,,yj)  ,^y(U,3(x+D^^^,y+d^^))A^^ 


i k 


exp{j  (2rf^^t+e^j^+(j)^j^(t)) } (40) 


where  as  in  Eq  (14),  mn  = N.  As  in  Eq  (22), 

- .^x(U.s(>t+Eik))  ^(U.s(y+dik>) 


where  the  shifting  theorem  of  Fourier  transforms  has  been 
used  (Ref  12:9;  27:65).  Now  Eq  (40)  can  be  written 
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Ui(xi,yi,t)  = P(xi  ,y,)Ujs(Xi  ,yj)  A.. 

i k 

exp{j  (2Tif^^tg+e^j^+<j)^j,  (t)) } (42) 


This  is  identical  to  Eq  (26) , the  field  for  the  multi-beam 

coaxial  case,  except  for  t where  t is  dependent  on  time 

s s 

and  space. 


t . t + £i!iL 


. + 

*ik  * ^Ik  * 

The  mean  of  this  field  is  similar  to  Eq  (27) 


(A3) 


E(U,(x^,yj,t))  = P(Xj.y,)U^g(x^.y^)  E E A.^ 


exp 


°ik 


2 '' 


(A4) 


and  the  covariance  is  similar  to  Eq  (28) 


Cov(Uj,U;]  = P(x^  ,yj)P*(Xj  .ypUjgCXj  ,yi)U*3(Xj  ,yp  E E A^^ 

X Ic 

exp  (i^Tif  i^(tg-tp)  {exp  (exp  [R())^^(t))  -1)  > (A5) 


where  t ^ is  a function  of  t' , x T,  and  y The  mean  of  this 

S ^ 1'  •'  1 

field's  intensity  is 


E(I.)  - 


TIW  Zi 


2Trb 

Xz 


r(xf+yj)] 

1 •' 


e(|£  E exp{j(2irf,^t^+e.^+*.^(c)))|=' 


(46) 


1 

t 

i 


i 


1 

i 

j 
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For  the  general  two  dimensional  array  this  expression  and  Eq 
(44)  cannot  be  evaluated  analytically  without  numerous  assump- 
tions, Eqs  (44^  and  (45)  are  similar  to  Eqs  (27)  and  (28) 


respectively  with  t replaced  by  t . Qualitatively,  the  two 
dimensional  array  produces  a pulse  train  in  the  far  field 
(there  are  still  N different  frequencies  being  added  so  the 
results  of  Eq  (31)  hold  for  = yj  = 0 in  Eq  (44)).  However, 
there  will  be  some  additional  spatial  modulation  because  of 
Eq  (43). 

Generalized  Results  and  Duality 

Eq  (42)  is  the  key  to  some  very  general  and  interesting 
results.  In  this  equation,  if  consid- 

ered independent  random  variables,  then  Eq  (44)  becomes 

E(  U j (Xj  ,yj  ,t)}  = P(Xj  ,yj)Ujs(Xj  ,yj)  E E A.y^  exp(je.,) 

i k 

(47) 

where  the  ^'s  represent  the  characteristic  functions  of  the 
indicated  random  variables  as  functions  of  the  indicated 
arguments  (i.e,  (a)  = E(exp(jax))  where  x is  a random 

variable)  and  is  evaluated  at  1.  Now  if  f,  D,  d,  and  (j) 
are  gaussian  (and  stationary)  with  means  f^^,  d^^  and 

and  variances  o^,  and  respectively,  then  Eq 

(47)  becomes  (see  Ref  26:159), 


E(Ui)  = P Ujg{exp  -%  C2iit)*o|+(2iT^^)®a^+(2ir^^)==o|+o^j } 


E 2 Aik  ^>‘PU<2''£ikts+®ik+*ik>) 


It  is  also  reasonable  that  = 0 (i.e.  f is  completely 
deterministic)  and  4)^^^  = 0 so  that 


E(Ui)  = P U,sexp{-% 


X 

Az 


(2it— i)2oJ+(27r-::-i-)=^a|+a 
1 


yiN2„2 .„2 


} Z E A k 
i k 


exp 


2tiD.,  X,  2TTd.,  y , 

j(2x£.^f+— + e..,.) 


Xzj 


Xz  j 


ik' 


(49) 


It  should  be  noted  that  Eq  (41)  is  still  valid  for  D and  d 
random  variables  because  the  Fourier  transform  is  taken  term 
by  term  and  each  Dj^k’  ‘^ik  nvmibers  (albeit  part  of 

a sample  space).  In  addition,  D,  d,  and  4)(t)  are  considered 
to  have  identical  variances  for  each  i,  k as  was  assumed  for 
Eq  (30). 


Now,  Eq  (49)  is  a useful  general  result.  If  i = k = 

Ojj  = = Dj^k  “ ~ ^ result  is  the  single  beam 

case  represented  by  Eq  (20)  with  = A and  6^^  =6.  If 
CTp  ~ ~ ^ik  ~ ^ik  ~ ^ result  is  the  coaxial  case  repre- 
sented by  Eq  (30)  with  A^k  “ 4.,  f^k  “ " ®j[  (i.e. 

the  double  series  collapses  to  a single  series) . If  = 

=0  the  result  is  the  general  two  dimensional  array  case 
with  Djj^k  deterministic  as  represented  by  Eq  (44) 

(with  o-k  = o^). 
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If  k = 0,  the  result  is  the  general  one  dimensional  linear 
array 


E(Ui (Xj ,yj ,t))  = P(Xj,y^)U  (Xj,y^)exp{-% 


) 


Z A.  exp 
i ^ 


2TfD.x, 

j(2iTf  t+  — -J— i-  +e.) 

^ AZi  ^ 


(50) 


The  deterministic  version  of  this  field  (a^  = 0)  will  be 
treated  in  the  last  section  as  well  as  the  case  for  D a 
random  variable. 

At  this  point,  three  souices  of  error  in  the  far  field 
that  are  a result  of  the  physical  positioning  of  the  laser 
elements  must  be  considered.  These  are  errors  in  the  element 
location  in  x (for  a linear  array,  x and  y for  a planar 
array) , errors  in  the  alignment  of  the  individual  optical 
axes  with  each  other,  and  errors  in  the  element  location  in 
z.  The  first  problem  can  be  considered  in  a manner  completely 
analagous  to  the  results  for  the  phase  variance  discussed 
in  connection  with  Eqs  (34) , (35) , and  (36) . For  the  linear 
array,  each  laser  may  have  an  error  represented  by  the  dif- 


ference between  its  actual  location  and  D. 


It  is  reasonable 


that  this  error  may  be  distributed  in  a gaussian  manner  so 
that  from  Eq  (50)  the  variance  of  that  error  is  o^.  Now  Ojj 
can  be  determined  as  a function  of  x,  and  z so  that  the 
result  has  a specified  effect  on  E(Uj),  For  example,  if 
if  Xj  = 1 m,  X = 10.6n,  z^  = 1 km,  and  ojj  = 1 mm,  then 


exp 


Xz, 


D 


= exp(-0.18)  =0.84.  For  a specified 
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amplitude  reduction  (for  example  807o  or  about  1 db)  and  a 
given  range  Zj,  the  accuracy,  o^,  with  which  the  laser  posi- 
tions must  be  m.ade  can  be  found.  The  result  is 

“d  - (51) 

where  ri  is  the  specified  reduction  and  a^,  Zj  , Xj , and  X are 
in  meters.  As  Zj  increases,  the  necessary  accuracy  of  the 
location  of  the  elements  decreases.  A similar  result  could 
be  derived  for  or  a planar  array  with  and  o^. 

As  for  Eq  (37),  the  covariance  of  the  field  varies 
with  element  location  as 


, 2ttx,  ■>  2itx. 

Cov[Uj  ,Uj)«exp|-(-^)^ay  (exp  | (-^)  2Rj^(x-x")  J -1)  (52) 


(Ref  26:159,476).  But  Rp(x-x')  = E(DD')=  E(d)e(D')  = Dj 
since  D and  D'  are  independent  (i.e.  the  location  of  one 
element  in  no  way  influences  the  location  of  the  next) . 
Since  the  element  locations  are  uncorrelated,  EfDD")  “ 
and  Eq  (52)  becomes 


cov(u,  ,u;)  cc 


1-exp 


. 2TrXi 
•(--^) 
Xzi 


(52a) 


as  -in  the  temporal  case  (although  Eq  (38)  is  a variance  and 
this  is  a covariance).  Thus  the  results  of  Eq  (49)  indicate 
the  complete  time-space  duality  of  this  problem.  The  statis- 
tical results  arc  similar  for  instabilities  in  time  or  space 
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(i.e.  positional  errors),  and  the  general  result  contains  all 
cases  of  interest. 

The  second  physical  positioning  problem  occurs  when  the 
optical  axes  of  the  lasers  are  not  parallel.  In  this  case 
Eq  (14)  becomes,  for  a linear  array, 

U,(x,y,t)  = E U (x+D.,y)  A 
i 

exp(j  (2Trfj^t+ej^+(f)j.  (t)  + ^ X sin  o^j^)}  (53) 

where  Oj^  is  the  angle  of  the  ith  optic  axis  with  respect  to 
the  normal  to  the  array.  For  small  angles,  sin  o^  o^  and 
the  far  field  becomes 

Ui(xj,yj,t)  = Uj  (Xj-2iTZjO^,yj , t)  (54) 

where  the  frequency  shifting  property  of  Fourier  transforms 
has  been  used  (Ref  27:65).  The  angle  o^^  now  affects  two 
factors:  exp  [ j y^(Xi -2itz  jO^)  from  Eq  (22)  and 
2TrD>  1 

exp  j — ^ (xi-27rz2  0.)  from  Eq  (42).  For  the  coaxial  case, 

I A Z J ^ j 

the  gaussian  envelope  behaves  as  in  Eq  (39),  therefore  vising, 
for  example,  the  Rayleigh  criterion  (Ref  15:354,360)  the 
maximum  o^  should  be  on  the  order  of  Wj  or  about  3 mrad  and 
no  larger.  Similar  results  are  obtained  for  y in  the  two 
dimensional  case. 

One  final,  rather  serious,  problem  in  the  area  of  laser 
position  errors  must  be  addressed.  Eqs  (26)  and  following 
were  derived  on  the  basis  that  all  N beams  that  were  super- 
imposed propagated  a distance  z^  into  the  far  field.  This 


J 
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assumes  that  Zj  is  measured  from  precisely  the  same  point  in 
each  output  beam  (for  example,  the  beam  waist).  If  this  is 
not  true,  i.e,  if  the  wa:ists  are  not  perfectly  aligned,  then 
the  Zj's  are  not  the  same  and  in  fact  must  be  represented 
as  Zi  +Az^  where  Az^  is  the  offset  of  the  ith  beam.  Substi- 
tuting this  new  distance  into  Eq  (26)  indicates  that  Az^ 
makes  a negligible  contribution  everywhere  except  the 
exp(jkAz^)  term  in  P(Xj,yj)  (see  Eq  (18)).  If  kAz^  is 
folded  into  the  interval  0 to  2ir  and  is  considered  a gaussian 
random  variable  independent  of,(i)^(t),  then  exp(jkAZj|^)  = 

f al) 

®xp  - ~ where  is  the  variance  of  kAz..  Even  if  the 

I 2 J z 1 

beam  waists  could  be  perfectly  aligned  so  that  Az^  was  zero, 
vibrations,  thermal  expansion,  and  other  effects  would  com- 
bine to  create  a Az^.  Presumably  this  additional  phase 
variation  could  be  combined  with  (|)j|^(t)  and  compensated  appro- 
priately. However,  even  this  requires  that  the  optical  path 
lengths  from  the  point  at  which  the  beams  are  sampled  to  the 
detectors  that  form  the  first  part  of  the  phase  locked  loop 
must  be  equal.  The  bottom  line  is  that  for  coherent  combin- 
ation, the  relative  phase  of  the  lasers  must  be  controlled 
precisely,  and  at  optical  wavelengths  that  is  a very  strin- 
gent requirement. 
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III.  Coaxial  Case 


The  coherent  combination  of  several  lasers  coaxially, 
i.e.  with  the  optical  axes  of  the  individual  lasers 
coincident,  led  to  a mean  far  field  given  by  Eq  (27). 
Certain  simplifying  assumptions  produced  Eq  (31) . This 
field  is  seen  to  be  a temporal  pulse  train  as  an  "envelope” 
on  a "carrier"  of  frequency  f - — ^ Af  where  f is  the 
original  optical  frequency.  In  addition,  the  signal  has  a 

f 2\ 

peak  amplitude  of  ANexp  - and  has  a spatial  distribu- 
tion  given  by  Eq  (23)  which  describes  a gaussian  spherical 
wave.  The  temporal  pulse  train  given  by 
characteristics  are  functions  of  N and  Af  suggests  that 
the  signal  could  be  used  for  ranging  (as  in  a radar) . If 
so,  a measure  is  needed  for  the  performance  of  such  a 
system  as  a function  of  controllable  parameters  such  as  N 
and  Af. 

Ambiguity  Function 

Before  any  performance  criteria  can  be  developed, 
there  must  be  some  assumptions  regarding  the  receiver 
structure.  Since  this  is  not  a receiver  design  study,  only 
general  assumptions  will  be  made.  (For  further  detail  see 
Ref  6,  11,  23,  32,  37,  and  40.)  The  first  assumption  of 
course  is  that  the  receiver  employs  coherent  (heterodyne) 
detection  (Ref  11:173-200).  Aside  from  the  requirement 
discussed  in  the  previous  section,  this  will  result  in 
better  performance  than  in  the  direct  detection  case  since 


both  amplitude  and  phase  information  are  preserved.  The 
heterodyne  detector,  however,  is  sensitive  to  spatial 
alignment,  but  this  is  not  assumed  to  be  a problem  here. 

The  second  assumption  is  that  a matched  filter  processor 
structure  is  used.  For  the  case  of  signals  in  additive 
white  gaussian  noise  (this  is  true  in  general  for  hetero- 
dyne receivers) , the  matched  filter  is  the  optimum  receiver 
(Ref  6:18-32;  23:335-395;  40:311-317).  The  matched  filter 
essentially  correlates  the  received  signal  with  a known 
signal  and  gives  a peak  response  when  the  two  signals  are 
"matched."  This  of  course  requires  that  the  receiver 
structure  change  if  the  signal  changes.  Fig.  5 is  a repre- 
sentative system  diagram.  The  third  assumption  is  that 
certain  physical  phenomena  that  affect  the  optical  signal 
as  it  is  transmitted  and  reflected  are  ignored.  These 


Combiner 


Local  Oscillator 
Field 


Fig.  5.  System  Diagram  of  Heterodyne  Receiver 
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include  such  effects  as  path  losses,  field  of  view,  losses 
in  the  optics,  field  alignment,  aperture  size,  and  others. 
These  are  ignored  because  they  affect  any  signal  in  the 
same  manner  and  therefore . are  unnecessary  when  comparing 
the  relative  performance  of  various  signal  structures . 

Based  on  detection  theory,  the  likelihood  ratio  is  a 
function  which,  when  compared  to  a threshold,  produces  an 
optimum  decision  of  which  signal  was  the  "most  likely" 
received.  This  likelihood  function,  in  the  case  of  a 
heterodyne  receiver  (and  gaussian  noise)  with  a matched 
filter  processor,  turns  out  to  be  proportional  to 

OO 

x(t,v)  = / s(t)  s*  (t-x)  exp  (j2Trvt)dt  (55) 

— CO 

which  is  just  a two  dimensional  correlation  of  a signal  s(t) 
(which  can  be  complex)  ii7ith  its  time  (x)  translated  and 
doppler  (v)  shifted  version  (Ref  6:59-108;  31:70-75,118-158; 
37:275-313).  The  ambiguity  function,  x(t,v),  is  a function 
of  time  delay,  x,  that  is  related  to  the  range  of  the  tar- 
get and  doppler  shift,  v,  that  is  related  to  the  velocity 
of  the  target.  (x  is  actually  "normalized"  so  that  x(0,0) 

corresponds  to  some  known  x and  v .)  The  signal  of 

0 0 

interest  here  is  just  the  mean  of  the  far  field  in  Eq  (31) . 
Since  the  only  part  of  the  signal  envelope  in  Eq  (31)  that 


AO 


varies  with  time  is  , and  since  the  detector  does 

sin  irAft 

the  spatial  integration  of  the  field,  the  signal  of  inter- 
est will  be  considered  to  be 


s(t)  = Aexp 


~T 


sin  NTrAft 
sin  irAft 


(56) 


where  Aexp 


is  retained  explicitly  as  part  of  the 


amplitude  of  s(t)  and  exp|^-  -g-j  is  considered  to  be  the 
strong  signal  approximation  to  the  characteristic  function 
of  the  partially  controlled  random  phase,  <(i^(t).  Unfor- 
tunately, substituting  Eq  (56)  into  Eq  (55)  yields  a 
result  that  is  not  immediately  integrable  to  obtain  x('f,v). 
However,  over  a single  period. 


0^1  sin  NirAft 

s(t)  = exp  - -J-  

\ J 


(57) 


This  approximation  gets  better  as  N gets  larger.  Therefore, 

,,\  _ A^  “ sin  NirAft  sin  NirAfCt-r)  c . o ..i  jx. 

- XF  -i  iFt FtFT ^ exp(j2iTVt)dt 

= ^ .^^(S(t)  S(t-T)) 

= ^ *^j-(s(t-T)) 

= ^ * exp(j2iTVT)J^^(s(t)}  (58) 
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Now,  by  the  duality  property  of  ambiguity  functions 
(Ref  37:309-310) 

X(t,v)  = Xj(v,-t) 

where  x is  the  ambiguity  function  of  s(t)  whose  Fourier 
transform  is  S^(f)  (whose  ambiguity  function  is 

Now  the  Fourier  transform  of  s(t)  in  Eq  (58)  is  just  a 
square  pulse  of  height  1 and  width  NAf  whose  ambiguity 
function  is  given  by  (Ref  37:280) 


IXj(t,v)1 


M] 


sin  NirAfv  1 


NA?J 


irv  1 


NAfJ 


, I T I 5NAf 


(59) 


1^0,  elsewhere 


Therefore, 


|x(t,v)1  = n^lXj(v,-T)| 


fi  - ^4] 

NAf 


sin  TrNAf(-T)[l  - M] 


lo. 


1 - -1-^ 
^ XT  * 


NAf 


Tr(-T)  ] 


sin  NirAfr 


NAl 


1 - 


ITT  1 


\ lv| 
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nai) 


, |v|<NAf 


(60) 


0,  elsev/here 
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Eq  (60)  is  sketched  in  Fig.  6 on  Page  43.  1x(t*0)I  is 

shown  in  Fig.  Ja,  and  IxCO.v)]  is  shown  in  Fig.  7b  found  on 
Page  43. 

It  is  now  possible  to  obtain  Cramer-Rao  lower  bounds 
for  estimates  of  t and  v in  terms  of  the  signal  energy  to 
noise  ratio,  the  bandwidth  of  the  signal,  and  the  duration 
of  the  signal.  These  are  estimates  of  delay  time  and 
doppler  velocity  for  a given  signal  structure  that  has  an 
ambiguity  function  x(''^>'^)*  The  (parameter)  estimate  of 
delay  is  denoted  t and  that  of  doppler  is  denoted  v.  The 
variances  (and  standard  deviations)  of  these  estimates  are 
a measure  of  the  precision  of  the  measurement,  i.e  a 
smaller  variance  implies  a more  precise  measurement.  Fol- 
lowing Ref  31:70-74  and  37:294-299  the  results  are 


\ “ (Var(T-T))^  > Y 

' 0 «' 


(61) 


(Var(v-v))^  t 


(62) 


where  is  the  standard  deviation  of  the  delay  estimator  t, 

is  the  standard  deviation  of  the  doppler  estimator  v, 

the  signal  energy  is  assumed  to  be  much  larger  than  the 

N 

noise,  the  noise  is  given  by  the  height,  of  its  power 

spectral  density  (the  noise  is  assumed  to  be  white  and 
gaussian),  and  y and  6 are 


Y*  = 4tt*  /“  fMS(f)12df  (63) 

6^  = 411*  /“  t*ls(t)|*dt 
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(64) 
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where  y represents  the  rms  bandwidth  of  the  signal  s(t)  and 
6 represents  the  rms  duration  of  the  signal.  (Eqs  (63)  and 
(64)  were  derived  in  Ref  31  and  37  for  normalized  signals 
s(t)  and  S(f).)  Eqs  (61)  and  (62)  also  assume  that  there 
is  no  coupling  between  y and  6,  i.e.  that  there  is  no  linear 
frequency  modulation  component  (Ref  31:74).  Now,  the  range 
precision  and  range  rate  precision  respectively  are  given 
by  (Ref  31:74) 


'R 


R 


c 

^ a 

2 V 


(65) 

(66) 


The  signal  energy,  E,  is  found  from  Eq  (56)  in  series  form. 


ls(t)  1 = exp|- 


sin  NirAft 


sin  irAft 


= exp 


0^'' 

X 


Z (exp  (j  2iTAf t)  ) ^ 
i 


Therefore , 

E = I s (t) I ^dt 
P 

where  T^  = ^ is  a period  (Eq  (33))  so  that 


] 


45 


E = 


1 

2Ar 

1 

■2Af 


Aexp 


Z (exp  (j  2'nAf  t) ) ^ 
i 


dt 


= A^exp 


= A^exp 


A^exp 


1 

211 

1 

2Af 


N + Z Z exp  ( j 271  (i-k)  Af  t) 


dt 


-0" 


-0^ 


fJL  4-  v V sin  Tr(i-k) 
(Af  7i(i-k)Af  _ 


N 

Af 


(67) 


Eq  (60)  was  derived  from  the  approximation  in  Eq  (57) . 
This  approximation  is  good  only  for  a single  period  of  the 
signal,  s(t)  (Eq  (56)),  However,  s(t)  is  periodic  in  t and 
can  be  represented  as 


m-1 

s(t)  = Z s^(t-iT^) 


i=0 


= ? ^a(t  - 
1 

where  m is  the  number  of  repetitions  of  s (t) , s (t)  is  the 

Si  d 

approximation  over  a single  period  , and  f p “ ^ • Now 


the  ambiguity  function  of  s(t)  becomes 

j7rv(m-l+p)^|x3[T-^,v] 


m-1 

x(t,v)  = Z exp 
p=-(m-l) 


sin  7rv(m-  Ip  I 


sin  TTV 


where  Xa(''^*^^  given  by  Eq  (60)  (Ref  31:185-186).  It  is 
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seen  that  the  new  ambiguity  function  is  just  the  sum  of  the 


i 


i 


old  ambiguity  functions  shifted  on  the  delay  axis  by 
intervals  of  ^ and  multiplied  by  a weighting  term.  The 
magnitude  of  the  pth  term  of  the  ambiguity  function  is 
just  (Ref  31:187), 


I 

1 : 


IxCt.v)  I 


at  p 


'a 


^ Af’'' 


sin  TTv(m-  jp  | 


sin  irv 


Af 


(68) 


Thus  there  are  delay  ambiguities  at  intervals  of  ^ and  a 

sin  Trv(m-  |p|)^ 

fine  structure  (due  to  j ) in  doppler  with 


sin  ttvt-c 
Af 

intervals  of  Af.  Now  some  specific  applications  for  a 
coaxial  system  whose  far  field  is  given  by  Eq  (31)  can  be 
considered. 


Range  Measurement 

When  only  range  measurements  are  desired,  the  doppler 
is  considered  to  be  known,  i.e.  tracked  (or  zero  for  a 
stationary  target).  If  the  doppler  is  unknown,  then  it 
must  be  considered  at  the  same  time  x is  considered  i.e. 
jointly  estimated  (this  is  done  later  in  this  section) , or 
the  ambiguity  function  must  be  broad  along  the  v axis  so 
that  there  is  little  chance  of  error  in  estimating  x. 
Therefore  the  interest  in  xC'^i'w)  is  just  along  the  x axis. 
Fig.  7a  is  the  x axis  for  a single  pulse.  Multiple  pulses 
(Eq  (68))  would  produce  a periodic  function  with  amplitude 
weighting.  As  long  as  the  round  trip  time  (delay)  of  a 
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pulse  is  less  than  the  pulse  rate,  there  is  no  ambiguity 
as  to  which  peak  on  the  t axis  is  involved.  The  width  of 
the  central  peak  in  Fig.  7a  is  a crude  measure  of  the 
precision  with  which  the  target  range  can  be  measured. 
Basically,  this  precision  is  a function  of  so  that  as 
the  total  spectral  width  of  the  field  is  increased  the 
range  precision  can  be  made  arbitrarily  good  within  the 
limit  that  NAf  is  less  than  or  equal  to  the  gain  line 
width  of  the  lasers.  Explicitly,  the  range  precision  is 
found  by  combining  Eqs  (61)  and  (65)  to  obtain 


Now  from  Eq  (63) , 


y"  = Att"  C fMS(f)rdf 


if  S(f)  is  normalized.  If  S(f)  is  not  normalized  then 


r fMs(f)rdf 


Now,  from  Eq  (57) 


_ * a^lsin  NirAft 

s(t)  = Aexp  -.J-  


Therefore , 


1 ' 


S(f) 
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Af 


Aexp  j^-- y 
[^0 , elsewhere 


< NAf 
“ 2 


= 4it" 


NAf 

• ^ 

1 

2 

2/  2 

0 

Af 

■ 

df 

df 


N^Af^ 


(71) 


So  that  using  Eqs  (67)  and  (71)  in  (69), 


12  exp  1 

1-cjMN 

N AJ 
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T 


(NAf)' 


C/3N7 


/Sir exp  (-a^  J N^Af 


(72) 


Eq  (72)  shows  that  the  range  measurement  precision  is  pro- 
portional to  which  is  in  line  with  what  was  expected. 


Since  NAf  is  the  total  bandwidth  of  the  signal  without 
regard  to  the  frequency  spacings,  the  assumption  of  f^  = iAf 
was  not  really  restrictive.  Eq  (72)  then  is  a rather  general 
result  that  represents  the  best  possible  range  measurement 
precision  obtainable  by  the  temporal  signal  in  the  mean  far 
field  of  Eq  (27). 


Velocity  Measurement 

In  the  dual  case,  a similar  derivation  can  be  made  for 
doppler  (velocity)  measurement  precision.  In  this  case  the 
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range  is  presumed  known  (or  the  ambiguity  function  is  broad 
in  t) , so  that  the  interest  in  x(t,v)  is  only  along  the  v 
axis.  Fig.  7b  is  the  v axis  for  a single  pulse.  Multiple 
pulses  introduce  a fine  structure  as  given  in  Eq  (68) . For 
a single  pulse,  the  crude  doppler  measurement  precision 


improves  as  NAf  gets  small.  This  implies  that  increasing 
2 

tp  = also  increases  the  doppler  precision.  Basically, 
this  says  that  the  more  information  that  is  received  (i.e. 


the  longer  the  pulse) , the  better  the  doppler  precision 
can  be  made.  From  Eqs  (62)  and  (66) 


R 


o /2l 


(73) 


and  from  Eq  (64), 


52  = 4tr^  /“  t^ls(t)|^dt 


so  that  for  the  unnormalized  s(t)  of  Eq  (57) 


6^  = C tns(t)|Mt 
/:is(t)rdt 


(74) 


00  M 

Therefore  for  a single  pulse,  /^|s(t)|^dt  = A*exp(-a^)^  as 


in  Eq  (67)  and 

1 

IZ  tMs(t)rdt  = 2 l^sln'NTTift  dt 


A^exp (-0^1 


f 1 


n^Af  ^ 


2Af 


(75) 


1 


50 


r 


vhere  the  integral  has  been  taken  over  a single  period  of 
the  signal.  For  m pulses,  the  signal  energy  is  multiplied 
by  m so  that  Eq  (74)  becomes 


6" 


Att' 


A^exp  f-o^l 


m 

24? 


fmA' 


exp^-o^')Nj' 


2 


(76) 


and  6 is  seen  to  be  independent  of  the  number  of  pulses  not- 
withstanding the  fine  structure  on  the  doppler  axis  given 
in  Eq  (68).  This  just  indicates  that  the  doppler  precision 
is  affected  only  by  the  pulse  width  and  not  the  number  of 
pulses.  The  fine  structure  in  Eq  (68)  tends  to  improve  the 
resolution  of  targets  with  similar  doppler  velocities,  but 
it  also  increases  the  number  of  doppler  ambiguities  (peaks) 
(Ref  31:191).  Eq  (76),  however,  is  concerned  only  with 
measurement  precision  and  not  resolution,  i.e.  local 
accuracy  without  regard  to  ambiguity.  Now  using  Eqs  (67) 
and  (76)  in  (73), 


i/Mi 

V N 


exp  f - o ] N 2 


Af  NAf2 


A/Nf.Af  ^ 

/16  A^exp{-o^) 


(77) 


and  the  doppler  precision  depends  only  on  Af. 

This  result  is  slightly  different  from  what  was  ex- 
pected from  inspection  of  Eq  (60)  since  N does  not  appear 


1 


I 
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in  Eq  (77).  The  reason  is  found  in  Eq  (75)  where  the  rms 
signal  duration  was  determined  for  a period.  When  this  is 
done  for  square  temporal  pulses  (as  in  radar  theory) , the 
result  is  an  integral  over  the  pulse  width.  If  this  is 
done  here,  by  analogy,  then  Eq  (76)  becomes 


6"  = 


(76a) 


so  that  Eq  (77)  becomes 


- 


XV'NoNAf  ^ 
/1 6 A^expf- 


(77a) 


The  discrepancy  occurs  because  s(t)  is  of  infinite  duration. 
This  is  the  dual  of  the  radar  problem  where  the  ideal  rec- 
tangular signals  have  infinite  mean  square  bandwidth 
(Ref  31:49-50).  In  that  case  the  signal  is  considered  band 
limited  and  appropriate  approximations  made.  The  result  is 
dependent  on  rms  pulse  width  as  it  is  in  Eq  (77a)  for  this 
case.  However,  for  the  lasers,  even  though  the  signal  is, 
strictly  speaking,  time  limited  (it  is  at  least  causal)  the 
infinite  result  may  be  more  closely  approximated  here  than 
for  the  radar  case.  Thus  Eqs  (77)  and  (77a)  represent  an 
interesting  dichotomy  of  results  that  are  both  reasonable, 
but  both  different.  These  equations  represent  the  best 
possible  precision  in  doppler  measurement  (when  range  is 
known)  that  can  be  obtained  from  the  field  of  Eq  (31) . 
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Coupled  Measurements 


It  is  often  desirable  to  make  measurements  of  both  range 
and  velocity,  i.e.  the  case  where  both  are  unknown.  Eqs  (72) 
and  (77)  or  (77a)  are  still  valid  as  long  as  there  is  no 
linear  frequency  modulation  on  the  signal,  i.e.  the  complex 
envelope  is  real  (Ref  31:74;  37:299).  If  this  is  the  case, 
the  range  precision  can  be  improved  by  increasing  N or  Af. 
However,  increasing  Af  degrades  doppler  precision  relatively 
rapidly.  By  increasing  N,  it  is  possible  to  improve  range 
precision  while  the  doppler  precision  remains  the  same 
(Eq  (77))  or  degrades  more  slowly  (Eq  (77a)).  Doppler  pre- 
cision can  be  increased  by  decreasing  Af,  but  to  maintain 
the  range  precision,  N must  be  increased  accordingly.  This 
is  a rather  nice  result  although  increasing  N becomes  pro- 
hibitively expensive  at  some  point  (even  for  the  Federal 
Government)  so  there  is  a limit  on  the  range  and  doppler 
precision  that  can  be  obtained.  In  addition,  since  NAf  must 
be  at  most  equal  to  the  laser  gain  line  width,  as  N is 
increased,  Af  is  decreased.  This  coupling  between  N and  Af 
further  complicates  the  task  of  obtaining  the  best  possible 
simultaneous  range  and  velocity  precisions. 


It  is  possible  to  improve  and  simultaneously  by 


E 


increasing  the  energy  to  noise  ratio,  If  is  fixed. 


this  can  be  accomplished  by  increasing  E.  From  Eq  (67),  E 
can  be  increased  by  increasing  N or  decreasing  Af  with  the 


same  effects  on  and  noted  above.  E can  also  be 


increased  by  decreasing  the  variance  of  the  phase,  i.e. 
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controlling  the  phase  more  tightly,  or  by  increasing  A. 

Both  of  these  approaches  can  be  taken,  but  both  have  limits; 
for  example,  a can  be  made  no  smaller  than  zero,  and  A is 
limited  by  the  size  of  the  laser  gain  medium  as  well  as 
other  characteristics  of  the  laser. 

The  simplest  way  to  increase  E without  changing  N , N, 


or  Af  is  just  to  collect  more  pulses.  Eq  (67)  was  developed 


for  the  energy  of  a single  period,  If  m pulses  are 


used,  the  energy  in  Eq  (67)  increases  linearly  with  m,  and 


and  from  Eqs  (72)  and  (77)  are  improved  accordingly. 


This  approach  is  only  limited  by  the  number  of  pulses  that 
can  be  received  and  by  the  assumption  that  the  target  range 
and  doppler  are  not  changing  appreciably  during  the  obser- 
vation period  (Ref  31:70). 

One  additional  case  of  interest  is  that  of  introducing 
some  sort  of  frequency  or  phase  modulation  on  the  signal 
s(t).  This  could  be  implemented,  conceptually  at  least,  by 
modulating  all  of  the  lasers  simultaneously  and  in  the  same 
manner.  A simple  case  is  that  of  linear  frequency  modula- 
tion where 


Sjjj(t)  = s(t)exp(j7rgt2) 


(78) 


By  a property  of  ambiguity  functions  (Ref  31:123;  37:290) 
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0,  elsewhere 


This  transformation,  in  effect,  shears  the  old  v axis  to 
V = gT . This  will  decrease  the  size  of  the  central  peaks 
on  the  T and  v axes  so  that  if  either  range  or  doppler  is 
known,  the  other  may  be  determined  with  better  precision 
than  before.  However,  for  the  coupled  case,  the  Cramer-Rao 
lower  bounds  are  increased  as  shown  in  Ref  31:73-74  and 
37:298,  Hence  a simple  linear  modulation  is  not  helpful  to 
the  coupled  measurement  case.  Much  more  complex  modulations 
such  as  pulse  coding  (Ref  37:314-323)  may  be  more  useful, 
but  their  implementation  for  N lasers  would  probably  be 
very  difficult. 
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IV.  Linear  Array 


The  mean  far  field  for  a planar  array  of  lasers  was 
given  by  Eq  (44) . This  equation  cannot  be  reduced  to  a 
convenient  closed  form  for  the  general  case.  It  can  be 
solved  numerically  for  a given  set  of 

’^ik’  some  simplifying  assumptions  are  needed  before 

any  further  analysis  is  possible.  The  first  such  assumption 
that  will  be  used  throughout  this  section  is  that  the  array 
will  be  considered  one  dimensional,  i.e.  a linear  array. 

Most  of  the  results  can  be  extended  to  two  or  more  dimen- 
sions if  desired,  and,  in  fact,  one  special  case  of  a 
planar  array  will  be  considered.  Nevertheless,  the  impor- 
tant results  can  be  demonstrated  with  a linear  array.  The 
second  assumption  used  in  this  section  is  that  the  locations 
of  the  lasers  are  precisely  known.  Therefore,  Eq  (50)  can 
be  rewritten  for  = 0 as , 


E(U  (x  y ,t))  = P(x  ,y  )U  (x  y )exp 

111  llAOll 


~T 


A^exp 


j[2TTf.t  + + e.' 

i ^ ij 


(80) 


2itD  . X 

Because  of  the  — . ^ ^ term  in  the  exponent,  this  mean  field 

A Z 

1 

has  a spatial  modulation  within  the  overall  gaussian  envelope 
described  by  P (x^ ,y ^ )U^ ^ (x^ ,y ^ ) . The  performance  measures 
in  this  section  will  deal  with  the  characteristics  of  this 
spatial  modulation  in  terms  of  N,  f^,  D^,  and  0^^^  and,  for 
certain  applications,  will  consider  peak-to-side  lobe  ratios. 
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Pattern  Characteristics 

In  order  to  maintain  the  pulsed  nature  of  the  signal 
used  in  the  coaxial  case,  the  first  assumption  used  here  is 
that  f^  = iAf;  other  cases  will  be  discussed  later.  In 
addition,  0^  is  a simple  phase  shift  that  merely  determines 
an  absolute  reference  for  either  the  time  or  space  coordi- 
nates. Hence,  with  no  loss  in  generality,  it  can  be 
assumed  to  be  zero.  An  interesting  use  for  9^  will  be 
discussed  later. 

Scanning . For  uniform  aperture  excitations  (A^  = A) 
and  considering  a uniform  array  where  = iD,  Eq  (80) 
becomes 


As  before,  the  exponential  factor  indicates  a shift  of  the 

optical  frequency  to  Af,  but  it  also  indicates  that 

there  is  now  a spatial  frequency  because  of  the  term. 

A z j 

The  ratio  of  sinusoids  in  Eq  (81)  can  be  written  as  a 

function  of  S(k"Xi+ut)  where  =.,-1^1  and  a)=27rAf.  This 
^ A J 

is  the  standard  form  of  a classical  travelling  wave  whose 

profile  is  given  at  t = 0,  temporal  period  is  , phase 

u 

velocity  is  , and  spatial  period  is  Therefore, 
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Temporal  period  “ T “ ^ ^ sec  (Eq  (33)) 


Temporal  pulse  width  = t = sec  (Eq  (32)) 


Af  Az 


Phase  velocity  ~ ^ “ — g~ 


Spatial  period  = = 


Spatial  pulse  width  = m J 

The  spatial  beam  profile  at  t = 0 is  identical  (within  a 

scaling  constant)  to  the  temporal  profile  at  x^  = 0 as 

shown  in  Fig,  4.  . Therefore,  in  addition  to  the  temporal 

pulse  train,  there  is  now  a spatial  pulse  train  scanning 

linearly  in  x^  at  a rate  of  = Af  Hz.  This  still  occurs 

P 

beneath  the  gaussian  envelope.  From  Eq  (39),  for  X = 10. 6y, 
Zj  = 1 km,  and  w = 1 mm,  the  waist  at  z^  is  about  3.5m. 

The  beam  width  is  equivalent  to  Xp  which  for  an  array  dimen- 
sion ND  = 1 m is  about  2 cm.  For  D = 1 cm  (this  would  be 
for  an  array  of  100  lasers  spaced  1 cm  apart) , the  spatial 
period  is  about  1 m.  The  result  is  illustrated  in  Fig,  8. 
The  small  beams  measured  by  Xp  will  move  right  or  left 
(depending  on  whether  f^  = iAf  or  (N  - i)Af,  i.e.  whether 
the  frequency  offsets  step  "up"  the  array  or  '.'down"  the 


array)  at  a speed  of 


and  the  pattern  will  repeat 


every  ^ sec.  Unfortunately,  Xp  places  more  stringent 
restrictions  on  the  misalignment  angles  of  the  laser  optical 
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axes  as  shovm  by  Eq  (54) . Now  the  maximum  should  be  on 

X 

the  order  of  or  about  10“ ® rad.  This  is  two  orders  of 
magnitude  more  restrictive  than  the  10“^  rad  for  the  coaxial 


case. 

This  is  a useful  result  for  a system  that  does  not 
require  resolution  of  close  targets  in  azimuth  to  less  than 
about  Wj  and  does  require  multiple  target  returns  (for 
example  to  reduce  glint  or  other  transient  phenomena). 
Unfortunately,  it  also  complicates  the  ambiguity  problem 
as  there  is  no  way  to  distinguish  which  beam,  x^,  is  pro- 
viding the  return.  Therefore,  the  resolution  capability 
inherent  in  the  small  Xp  is  being  wasted.  These  side 
lobes  away  from  the  x,  = 0 peak  are  called  "grating  lobes" 


N Af  Xz 


1 m 


\ Gaussian 
\ Envelope 


i !i  I!  Offl  l ' 


Fig.  8.  Example  of  Far  Field  Mean  Beam  Pattern 
for  N=100,  D=1  cm,  Zj=l  km. 
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in  radar  theory.  It  can  be  sho^-m  that  these  lobes  do  not 
exist  if  D<^  (Ref  35:90)  which  is  impossible  at  optical 
wavelengths.  Eqs  (39)  and  (82)  suggest  two  possible  solu- 
tions. If  the  waist,  w^,  could  be  made  sufficiently  small 
(<Xp)  then  the  gaussian  envelope  could  be  made  to  attenuate 
unwanted  lobes.  From  Eq  (39)  the  only  way  to  do  this  is  to 
increase  w.  For  significant  attenuation  of  the  first  side 
lobes,  w,  - ^ , therefore  w ^ ~ . This  is  not  desirable 
since  originally  the  laser  aperture  radius  was  chosen  to 
be  to  minimize  diffraction  effects  and  maximize  power 
transmission.  Hence  the  aperture  radius  must  be  at  least 
on  the  order  of  the  aperture  separation,  and  the  apertures 
overlap.  Eq  (82)  indicates  that  can  be  increased  by 
decreasing  D and  similar  results  are  obtained,  i.e.  the 
apertures  overlap.  Thus  for  D>^,  there  is  only  one  other 
way  to  eliminate  the  grating  lobes  and  that  is  to  use  an 
aperiodic  array,  i.e.  D^^iD.  This  case  will  be  discussed 
at  the  end  of  this  section. 

Pattern  Shape.  So  far  the  amplitudes,  A^,  of  the 

lasers  have  been  assumed  to  be  equal  and  the  result  has 
S XTl  Nx 

Deen  a ■ ■■■. . ^ pattern  in  space,  Eq  (81).  It  is  useful  to 

pursue  a slightly  different  analysis  of  Eq  (80)  for  two 
reasons:  1)  more  flexibility  may  be  obtained  if  the  are 
not  all  the  same,  and  2)  the  side  lobe  structure  of  Eq  (81) 
may  not  be  desirable.  In  Eq  (81),  the  first  side  lobe  is 
the  highest.  Fig.  9 shows  the  peak- to-first-side  lobe 
ratio  as  a function  of  N.  Therefore,  even  as  the  number 


60 


N 


Fig.  9.  Peak- to-First-Side  Lobe  Ratio  versus 
N for  sin(Nx) /sin(x) . 


of  lasers  increases,  after  N » 20  the  highest  side  lobe 
grows  linearly  with  the  peak. 

If  Eq  (80)  is  rewritten  for  f^  = iAf,  = iD,  and 
0^  = i0  the  result  is 

E(U  (x  ,y  ,t))  = P(x  y )U  (x  ,y  )exp[-^]zA.z^  (83) 

111  lll^ll  ^ ^ X 


where 


z = exp 


211  Dx 

j (2irAft  + ■ + 0) 


(84) 


and  the  series  in  Eq  (83)  can  be  written 


lA. = A.,  , n (z-z.) 
^ i N-1  X 


(85) 


which  is  just  a polynomial  in  z (Ref  35:87).  Eq  (84)  is 
referred  to  in  antenna  theory  as  a z transform.  (For  a 
more  general  discussion  of  antenna  theory,  see,  for  example, 
Ref  5;  properties  of  z transforms  and  their  uses  may  be 
found  in  Ref  25:45-77  and  35:87-111.)  Now,  Eq  (84)  can  be 
written 


j 


1 
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r 

i 


z 


exp 


exp 


. 

« 

- sin  ej 


(86) 


X 

where  ^ in  the  far  field  approximation  is  the  angle  of 
observation  with  respect  to  the  z axis  (optical  axis)  and 
represents  the  angle  at  which  the  "main  lobe"  is  point- 
ing (Ref  35:89)  as  shown  in  Fig.  10. 


Of  course,  |sin  0^1  cannot  be  greater  than  1 so  that  values 
of  that  are  greater  than  1 essentially  repre- 

sent beam  angles  outside  of  the  "visible"  region.  In  fact, 
|X^t  ^ may  be  considered  modulo  2tt , i.e. 

2mr  + sin  0p  where  n = 0,  ±1,  ±2,...  and  |sin  0^1  < 1 
(Ref  35:89).  Therefore,  in  Eq  (86)  the  main  beam  location 
is  a function  of  time  as  was  shown  by  Eq  (82) . It  is  now 
possible  to  choose  the  zeroes  and  amplitude  weightings  of  the 
polynomial  in  Eq  (85)  to  yield  desirable  beam  width  or 
side  lobe  properties  (see,  for  example.  Ref  30,  36,  and 
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Ref  35:95-110).  The  array  can,  in  fact,  be  designed  so 
that  the  side  lobes  are  always  below,  a given  design  value 
(which  can  be  made  much  lower  than  the  results  in  Fig.  9) 
for  a given  main  beam  width  (this  is  called  a Chebyshev 
array;  Ref  35:111-118).  Thus,  with  this  approach  the  side 
lobes  may  be  better  controlled,  however,  there  are  still 
grating  lobes  because  b>2. 

It  is  interesting  to  note  that  by  duality  similar 
results  can  be  obtained  for  the  temporal  pulse  train  rep- 
resented by  Eq  (31).  In  this  case,  for  = 0,  the  zeroes 
of  Eq  (85)  may  be  chosen  so  that  the  temporal  side  lobes 
are  controlled  in  a desired  manner.  Presumably  appropriate 
z^  and  could  be  chosen  to  control  the  temporal  and 
spatial  side  lobes  simultaneously. 

Monochromatic  Beam 

The  z transforms  also  indicate  another  interesting 
possibility.  If  in  Eq  (86),  Af  = 0,  then  the  beam  position 
is  not  a function  of  time;  it  is  only  a function  of  9. 

This  represents  the  case  of  a continuous  wave  resultant 
beam  where  all  N lasers  are  oscillating  at  the  same  fre- 
quency, and  are  being  added  coherently  with  a control- 

lable phase  offset,  0,  across  the  array.  This  case  is 
very  similar  to  the  phased  array  in  antenna  theory. 

Beam  Steering . For  Af  = 0,  Eq  (86)  becomes 

z = exp  - 0o]  (87) 
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X 0 

where  “Qjj”  since  0^  is  small.  For  the  linear  array. 


then,  the  main  beam  may  be  pointed  anywhere  within  the 
gaussian  envelope  by  proper  choice  of  0.  (Of  course, 
because  of  the  grating  lobes,  one  "main  beam"  is  indis- 
tinguishable from  its  neighbor  and  in  fact  all  of  the  lobes 
move  simultaneously  with  0.  Nevertheless  a "main  beam" 
is  referred  to  for  convenience.)  If  desired,  the  main  beam 
may  be  made  to  scan  sinusoidally  for 

n - X0 

0 2itD 

= rsin  2TTf  t 
s 

0 = sin  2Trfgt  (88) 


where  r is  a scale  factor  that  determines  the  limits  of  the 
scan  displacement  and  f^  is  the  sinusoidal  scan  rate.  Two 
such  arrays  arranged  orthogonally  would  provide  two  sets  of 
pencil  beams  that  scan  orthogonally.  Such  a system  could 
be  used  for  tracking  a target. 

Planar  Array . For  the  special  case  of  f^  = Af  = 0, 

^ik  °ik  ^ ^ik  ^ ‘^ik 

written  as 


exp(j(^)(D^d^)+Je.J 


Now,  if  0^j^ 


ia  + ka  , Eq 

X y 


(89)  becomes  "separable". 
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expf  ji(^)D( 


V Xot  ' 
z 2irD'_ 


exp(jk(^)d(|l<-2^)] 


In  spherical  coordinates,  - z^sintjjcos^  and  = 
z^sinV^sin?  where  ip  is  the  angle  off  of  the  normal  to  the 
array  and  c is  the  angle  in  the  plane  of  the  array  (Fig.  11) 
The  result  is  the  same  as  in  Ref  10:36  (see  also  Ref 
5:141-142,  172-173).  Therefore,  for  a symmetric  distribu- 
tion of  lasers  around  the  center  of  the  array  and  with  a 


Direction  of 
Main  Beam 


• = Lasers 


Fig.  11.  Geometry  of  Planar  Array  and  Far  Field 

laser  at  the  center  of  the  array  (i.e.  an  odd  nximber  of 
lasers  in  both  rows  and  columns  which  would  require  at 


i 

J 


! 

i 


least  9 lasers)  the  pointing  angle  of  the  main  beam  is  at 


L 


ij;  , where  (Ref  10:38) 
0 0 


tan  c 

0 


D 


sin^  \1)  = 

0 


(91) 

(92) 


A conical  scan  is  useful  for  tracking  applications.  In 

terms  of  4*  • C a conical  scan  would  require  a constant  ip 
0 0 0 

and  a C that  increases  linearly  with  time.  If  D = d, 

0. 

then 


a = rcos  2Trf„t 
X s 

(93) 

tty  = rsin 

will  produce  a conical  scan  at  a rate  of  f ^ , For  this 

case,  4^^  = sin  * and  = 2irfgt.  Because  of  the  grating 

lobes,  the  resulting  scan  in  the  x , y plane  will  look  like 

1 1 

a series  of  concentric  circles  with  amplitudes  determined  by 
the  overall  gaussian  envelope. 

Aperiodic  Arrays 

As  noted  earlier,  the  grating  lobes  may  not  be  objec- 
tionable depending  on  the  application.  However,  they  must 

be  removed  to  obtain  the  resolution  inherent  in  x and  to 

P 

eliminate  ambiguities  in  target  location  due  to  the  numerous 
lobes.  The  most  feasible  way  to  remove  them  is  to  use  an 
aperiodic  array  where  D^^^iD.  The  analysis  here  is  similar 
to  radar  antenna  theory  (Ref  5,  30,  35)  except  that  the 
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wavelengths  are  much  shorter.  In  general,  aperiodic  arrays 
can  be  designed  with  the  element  spacings  chosen  according 
to  an  algorithm  or  chosen  in  a completely  random  manner. 

The  statistics  for  an  ensemble  of  random  array  configura- 
tions lead  to  various  descriptions  of  the  far  field  beam 
pattern.  Both  types  of  aperiodic  arrays  are  discussed 
briefly  with  general  results  dra\TO  from  antenna  theory  and 
applied  to  the  laser  case. 

Spatial  Taper.  Unfortunately  there  is  no  general  theory 
governing  the  design  of  aperiodic  arrays.  Several  techniques 
are  listed  in  Ref  35:135;  the  most  successful  of  which  is 
apparently  a process  called  dynamic  programming  which  is  a 
sophisticated  trial-and-error  procedure.  Iterative  methods 
have  also  been  used  to  synthesize  a desired  pattern  (Ref  36) . 
Ref  20  showed  that  for  a linear  array,  the  lowest  side  lobe 
levels  were  generally  from  space  tapered  arrays.  Spatial 
tapering  usually  results  in  an  array  in  which  the  radiating 
elements  are  closely  spaced  near  the  center  of  the  array  and 
become  more  spread  out  near  the  edge  of  the  array.  For 
example,  the  element  spacings  can  be  made  inversely  propor- 
tional to  the  aperture  excitation  (which  is  considered 
continuous  over  the  entire  array  aperture)  (Ref  35:126-129). 
The  result  is  a radiation  pattern  whose  main  beam  and  first 
couple  of  side  lobes  approximate  the  design  goals  fairly 
closely.  However,  the  side  lobes  tend  to  get  higher  for 
angles  further  from  the  main  beam  (Ref  35:130).  Fig.  12  on 
the  next  page  is  a representation  of  the  radiation  pattern 


67 


Fig.  12.  Representation  of  Typical  Space  Tapered 
Array  Pattern  (Adapted  from  Ref  35:130). 

for  a typical  spatially  tapered  array.  Region  I approxi- 
mates the  design  goals  with  low  side  lobes.  Region  II  is 
a transition  region  where  the  side  lobes  get  higher. 

Region  III  has  average  side  lobe  levels  similar  to  the 
results  for  random  arrays.  If  the  pattern  in  Regions  II 
and  III  will  not  be  used,  then  perhaps  spatial  tapering  can 
eliminate  the  grating  lobe  problem.  In  the  case  of  optical 
frequencies,  this  might  be  true  because  the  overall  gaussian 
envelope  severely  attenuates  the  pattern  away  from  the  main 

beam.  For  example,  at  a distance  of  4w  from  x =0,  the 

11 

gaussian  profile  is  about  36  db  do^^n  from  the  peak  at  x^  = 0. 
Unfortunately,  because  of  the  small  wavelength  involved. 
Region  III  occurs  much  closer  to  the  main  beam  than  for  the 
radio  frequency  case.  In  fact,  for  most  spatial  tapers,  the 
design  pattern  approaches  the  Region  III  results  by  the  first 
side  lobe  (Ref  35:150).  Since  the  location  of  the  first  side 
lobe  is  proportional  to  ^ where  L is  the  overall  length  of 
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the  array,  the  gaussian  envelope  will  not  be  selective 
enough  in  the  neighborhood  of  the  main  beam  to  kill  off  the 
high  side  lobes  in  Regions  II  and  III. 

Random  Arrays . Strictly  speaking  there  is  no  such 
thing  as  a random  array  since  any  physical  array  is 
completely  deterministic  (not  including  measurement  errors) . 
This  seems  to  imply  that  to  determine  anything  about  the 
far  field  pattern  of  an  infinitely  large  class  of  specific 
array  configurations,  one  must  realize  each  array  and  com- 
pute its  pattern.  Fortunately  this  is  not  the  case.  It  is 
possible  to  make  certain  probabilistic  statements  that  are 
valid  for  such  a class  of  arrays  whose  elements  may  be 
positioned  in  some  random  manner.  Such  statements,  though, 
are  merely  existence  statements  or  bounds  in  that,  within 
the  class  of  arrays  being  considered,  at  least  some  array 
configurations  exist  whose  patterns  are  described  by  these 
probabilistic  .statements.  This  approach  is  absolutely  no 
help  at  all  ir,  realizing  a specific  array  with  the  desired 
characteristics.  Such  a realization  must  be  left  to  a 
computer  or  a designer's  insight.  Nevertheless,  it  is 
useful  to  review  the  results  of  random  array  analysis 
(Ref  19,  20,  30,  35). 

Eq  (50)  was  the  result  for  the  general  one  dimensional 
’ ir  irray  wb*-rc  4^(t)  and  D were  considered  gaussian 

.*  ' ' When  considering  the  far  field  pattern, 

, t ir  .irbitr.iry  and  will  be  set  to 


E(U/x^,y^.t))  = P(x^.y^)U^g(x^.y^)exp  -%[0.^2+[2Tr^-2^j  a^^jj  E 


2TrD.x 


where  the  are  the  means  of  the  element  locations , and 
Ojj*  is  their  variance  (since  they  are  considered  to  have 
identical  variances).  If  the  array  is  uniform  (D^  = iD) , 
the  series  in  Eq  (94)  has  the  spatial  grating  lobes  shown 
in  Fig.  8.  It  can  be  seen  imaiediately  that  for  x^  0, 
can  be  appropriately  chosen  to  attenuate  the  grating  lobes. 
In  Eq  (51)  , v;as  chosen  small  so  that  the  far  field 
pattern  was  not  significantly  affected;  here,  however, 

(X  ' ^ 

relatively  large  to  accomplish  the 

desired  result.  In  a manner  completely  analagous  to  the 
temporal  results  of  Eqs  (34)  and  (35)  , the  mean  far  field 
intensity  is  proportional  to  N for  x^  # 0 since 

-v  sin  ^ 

E(1  J.  l-exp|2„^T  N + exp|2T^?  a "]  (95 

o-  !■'  JJ  I'-  I-*  J Sin  TTDr-zi- 


where  = A and  = iD  in  Eq  (94).  At  x^  = 0,  however, 

the  main  beam  peak  is  proportional  to  . Therefore,  the 

power  ratio  of  the  average  side  lobe  level  to  the  main  beam 
N 1 

peak  is  jjr  = jj.  so  that  the  average  side  lobe  level  is 
inversely  proportional  to  the  number  of  lasers  in  the  array 
(Ref  35:140-142). 
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There  is  still,  however,  a finite  probability  that 
some  side  lobe  will  be  much  higher  than  the  average  side 
lobe  level  of  This  is  true  since  the  far  field  pattern 
of  any  given  array  is  simply  one  possible  sample  function 
of  all  of  those  under  consideration,  and  any  individual 
sample  function  may  have  very  high  or  very  low  side  lobes. 

Eq  (94)  indicates  that  the  mean  of  can  be  made  to  atten- 
uate the  grating  lobes  for  x^/0  as  much  as  desired.  The 
field,  however,  may  deviate  from  this  mean  by  a considerable 
amount.  , To  investigate  this  deviation  a nev?  definition  is 
used  so  that  the  field  under  consideration  is  zero  mean 
and  unit  variance 


U (x  ,y  ,t)  - E(U  (x  ,y  ,t)} 
U"(x  ,y  ,t)  = i 1 5 5 i i 
^ ^ ^ vA/ar(U^  ,U^“) 


(96) 


In  general,  U'  is  a complex  random  process  with  real  and 
imaginary  parts.  For  a large  number  of  elements,  these 
parts  may  be  considered  asymptotically  jointly  gaussian 
random  processes  by  the  central  limit  theorem  (Ref  19:258; 
35:147).  With  this  assumption,  the  joint  probability 
density  of  the  real  and  imaginary  parts  of  the  field  U'  is 
known,  and  the  probability  that  any  side  lobe  (a  given 
amplitude  of  the  field)  will  not  exceed  some  specified 
threshold  can  be  determined.  (This  is  not  the' same  as  the 
asstimption  that  was  made  earlier  that  only  the  statistics 
of  the  phase  were  known  and  that  the  amplitude  was  deter- 
ministic. In  this  case,  both  the  amplitude  and  phase  of 
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the  far  field  U"  are  considered  random  for  this  class  of 

1 

arrays  since  the  elements  may  have  locations  determined  by 
some  unspecified  probability  density.)  Since  is  zero 
mean,  its  real  and  imaginary  parts  are  identically  distri- 
buted as  well  as  independent.  This  is  because  the  element 
locations  are  considered  to  have  a common  underlying  proba- 
bility density,  so  that  the  element  positions  are  identically 
distributed  about  their  means  D^.  Now  Var(U^U^*]  = N from 
Eq  (95)  since  E(lJ  = E(U^U^*)  = Var[U^U^*)  and  is  big 
enough.  Following  Ref  19:261  and  35:147-151,  the  result  is 
as  follows : 


A ' 

A ' 

Pr 

= 1-Pr 

/NJ 

^ Pr 


j 


exp 


exp 


exp(-B) 


(97) 


where  B = is  the  power  ratio  of  the  desired  threshold 
to  the  average  side  lobe"  level.  The  exponential  dependence 
is,  obtained  from  the  exponential  distribution  of  the  ampli- 
tude of  lU^|.  The  probability  that  the  peak  side  lobe  is 
less  than  B is  just  1-e  . Eq  (97),  however,  is  true  for 

only  a single  sample  value  of  at  some  x . To  find  the 
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probability  that  no  side  lobe  anywhere  in  the  pattern 
exceeds  B,  it  is  necessary  to  take  many  spatial  samples  of 
U'.  The  issue  now  is  how  many  samples  should  there  be  and 
where  should  they  be.  The  latter  problem  is  reduced  by 
simply  choosing  the  distance  between  samples  to  be  the  same. 
To  make  the  computations  simpler,  the  number  of  samples  may 
be  constrained  by  the  requirement  that  the  samples  be 
statistically  independent.  Crudely,  then,  the  number  of 
samples  needed  to  describe  the  random  process  is  a set  of 
independent  random  variables.  For  such  a case,  the  samples 
are  separated  by  no  less  than  the  coherence  distance  of  the 
pattern  and  the  pattern  is  considered  piecewise  constant 
over  a coherence  length.  The  total  probability  that  the 
set  of  {lU'(Xj^)|}  is  less  than  B is  just  the  product  of  the 
probabilities  that  each  individual  sample  is  less  than  B, 

B = Pr({|U'(xp|2}<B) 

= Pr(|U^(x^) |2<B  and  |U^ (x^) | 2<B  and. . . |U^(x^)|*<B) 

= (1-exp  (-B))”  (98) 


where  n is  the  number  of  samples  and  is  called  the  array 
parameter.  As  the  coherence  distance  of  the  pattern  gets 
small,  n gets  large.  Eq  (98)  can  be  rewritten, 


B = 


In 


1-  3^1 


(98a) 


Now  if  the  probability  (confidence  level)  3 is  chosen  and 
n is  known,  then  the  peak  side  lobe  level  (with  respect  to 
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the  average  side  lobe  level,  can  be  determined.  For 

example,  for  a 90%  probability  that  no  side  lobe  exceeds  B 
and  for  n = 10,  then  B = 4.56,  i.e.  with  90%  confidence, 
the  peak  side  lobe  of  |U^|^  is  4.56  times  as  high  as  the 
average  side  lobe  level,  i.  It  should  be  noted,  too,  that 
the  underlying  distribution  of  the  array  element  locations 
is  not  a factor  in  this  analysis;  this  is  one  of  the 
reasons  why  the  analysis  is  general  for  a large  class  of 
array  configurations. 

The  determination  of  the  array  parameter,  n,  can  be 
made  in  several  ways  (Ref  35:151-156)  and  the  result  can  be 
given  by 

n = ^ (99) 

where  L is  the  length  of  the  array.  This  is  intuitively 

correct  since  ^ is  roughly  the  angular  lobe  width  of  the 

pattern,  therefore  the  necessary  "sample  rate"  should  be 

proportional  to  For  optical  wavelengths,  then, 

(particularly  for  CO  lasers)  n is  on  the  order  of  10®.  The 

2 

effects  on  n of  a main  beam  at  some  location  other  than 
x^  = 0,  unequal  element  excitations,  non- isotropic  elements, 
and  signal  bandwidth  are  all  negligible  (to  within  an  order 
of  magnitude)  for  this  case  (Ref  35:171-184). 

Eq  (98)  overlooked  a minor  problem.  Since  B itself 
has  a variance  and  the  samples  are  not  infinitely  close 
together  (continuous),  the  actual  peak  side  lobe  is  higher 
than  B with  probability  one,  i.e.  U'  may  not  be  piecewise 
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constant  over  a coherence  distance.  Ref  19  ignores  the 
error  and  Ref  35  presents  a rather  obtuse  derivation  of 

2 

an  additive  correction  factor  which  turns  out  to  be  1 + g, 

2 

so  that  the  actual  peak  side  lobe  is  given  by  B + 1 + ^. 
For  the  optical  wavelengths,  the  correction  factor  is  neg- 
ligible, since  n = 10®  and  Eq  (98a)  gives  B s 12  (for  a 
confidence  level  of  50%)  and  B = 14  (for  $ = 90%).  Thus, 
the  correction  affects  the  actual  peak  side  lobe  less  than 
about  10%  numerically,  and  in  fact  merely  represents  a 
small  decrease  in  the  confidence  level  6.  Since  B is  the 
ratio  of  the  peak  side  lobe  to  the  average  side  lobe,  the 
ratio  of  the  peak  side  lobe  to  the  main  beam  is  given  by 


N 


In 


1-3 


(100) 


The  prv'.mary  restriction  on  Eq  (100)  results  from  the 
original  assumption  that  the  field  had  real  and  imaginary 
parts  that  were  asymptotically  jointly  gaussian.  Ref  19:262 
notes  that  this  approximation  is  valid  for  <<1.  A plot 
of  N versus  h is  given  in  Fig.  13  for  q = 1,5  where  L = 10*^X 
and  3 = 90%.  The  restriction  — tt  <<  1 is  shown  by  a dotted 
line.  Left  of  this  line,  the  confidence  level  may  be  less 
than  90%.  The  figure  is  very  nearly  the  same  as  that  in 
Ref  19:262.  The  figure  indicates  that  as  L(q)  increases, 
the  peak  side  lobe  is  also  likely  to  increase,  but  not  a 
great  deal  relatively. 


I 
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Fig.  13  also  indicates  a very  striking  result.  For  a 
reasonable  peak  side  lobe  level  of  -20  db,  aim  array  (q-5) 
must  consist  of  about  1400  lasers.  As  noted  earlier,  these 
probability  statements  are  existence  statements  for  an 
infinitely  large  class  of  arrays  whose  elements  are  positioned 
according  to  some  underlying  probability  density.  Therefore, 
there  are  at  least  some  arrays  which  have  about  1400  lasers 
whose  side  lobes  are  less  than  -20  db  with  a confidence  level 
of  about  907o.  This  is  not  true  for  most  of  the  arrays  in  the 
• class.  Furthermore,  this  is  not  a design  procedure  for  a 


specific  array.  Finding  an  array  with  the  above  performance 
is  by  no  means  a trivial  task  (analytically,  it  is  totally 
intractable).  However,  even  disregarding  the  difficulty  of 


designing  the  appropriate  array,  there  is  the  not  inconsid 
erable  problem  of  actually  building  it. 


V.  Conclusion 


k 


Summary  of  Results 

In  the  first  section,  classical,  amplitude  stabilized 
fields  were  used  to  describe  the  laser  output.  The  fields 
were  considered  complex,  coherence  separable  random  processes. 
The  Huygens-Fresnel  principle  was  used  to  propagate  the  field, 
and  the  far  field  (Fraunhofer  condition)  was  considered.  The 
general  result  for  the  superposition  of  N fields  contained 
all  of  the  cases  discussed,  and  the  time-space  duality  of  the 
calculations  was  evident.  When  coherent  detection  was  used 
to  measure  a field  at  any  space- time  point,  the  resulting 
measured  phases  of  the  fields  were  considered  stationary, 
gaussian  random  processes.  The  effect  on  the  far  field  pat- 
tern was  an  amplitude  reduction  proportional  to  the  variance 
of  the  phase.  As  the  variance  increased,  the  result  approach- 
ed the  incoherent  case;  and  as  the  variance  approached  zero, 
the  fields  we"^®  ceen  to  add  coherently.  A phase  locked  loop 
was  postulated  to  lock  each  laser  to  a reference.  A local 
oscillator  could  be  used  in  a coherent  (heterodyne)  detector, 
and  a phase  locked  loop  could  lock  the  reference  to  the 
local  oscillator.  The  receiver  then  could  track  the  signal 
without  regard  to  the  absolute  phase  variations  of  the  sig- 
nal or  the  local  oscillator.  Three  cases  of  misalignment 
were  considered.  Variations  in  laser  spacing  produced 
results  similar  to  the  effect  of  the  phase  variance.  In 
each  case  the  required  standard  deviation  for  a desired 
effect  on  the  far  field  was  determined.  Misalignment  of  the 
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laser  axes  was  also  considered  and  found  to  require  alignment 
accuracy  to  within  milliradians  for  the  coaxial  case.  One 
other  problem  investigated  was  that  of  the  misalignment  of 
the  fields  in  the  direction  of  propagation.  This  was  seen 
to  produce  a non-negligible  phase  variation  in  the  far  field 
that  required  physical  alignment  of  the  lasers  to  within  an 
optical  wavelength  to  eliminate  the  variation.  The  sta- 
tistics resulted  in  another  amplitude  reduction  of  the  far 
field  proportional  to  the  variance  of  the  misalignment. 

A few,  minor,  simplifying  assumptions  in  the  first  sec- 
tion produced  a far  field  result  almost  identical  to  that  of 
a mode- locked  laser.  However,  there  was  additional  flexi- 
bility in  that  the  pulse  widths  and  pulse  rates  could  be 
easily  varied  by  changing  the  frequencies  of  the  individual 
lasers.  The  minimum  pulse  width  was  determined  to  be  about 
3 nsec  using  COj  waveguide  lasers.  A laser  array  also  pro- 
duced spatial  modulation.  In  both  cases  the  overall  beam 
propagated  as  a gaussian  spherical  wave. 

In  the  second  section,  the  case  of  several  lasers  whose 
optical  axes  were  coincident  was  considered  in  more  detail. 
For  the  reasonable  assumptions  of  coherent  detection  and  a 
matched  filter  receiver,  the  performance  of  the  mean  of  the 
far  field  was  seen  to  be  described  in  terms  of  the  ambiguity 
function.  The  general  characteristics  of  the  ambiguity 
function  depended  essentially  on  the  total  bandwidth  of  the 
combined  field.  Thus  using  other  than  a constant  frequency 
interval  Af  would  not  be  significantly  helpful.  The 
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measurement  precisions  for  range  and  velocity  were  found  to 
be  proportional  to  and  /l^f ^ respectively.  Therefore, 

in  theory,  the  measurement  precisions  could  almost  independ- 
ently be  made  as  small  as  desired  within  the  constraint  that 
NAf  must  be  less  than  or  equal  to  the  laser  gain  line  width. 
Multiple  pulses  could  also  be  used  to  increase  the  signal 
energy  to  noise  ratio  thus  improving  both  measurements. 

Overall  modulation  such  as  linear  frequency  modulation  (chirp) 
was  not  seen  to  be  particularly  useful  in  improving  measure- 
ment precision.  This  coaxial  case  does  have  application  to 
target  detection  and  tracking  problems  as  well’ as  terrain 
mapping.  The  system  is  very  flexible  because  of  the  capa- 
bility to  change  pulse  characteristics,  thus  changing  meas- 
urement precision  when  necessary. 

In  the  third  section,  the  case  of  an  array  of  lasers 
was  considered.  A linear  array  was  analyzed  and  found  to 
have  a spatial  pulse  that  scanned  linearly  beneath  the  over- 
all gaussian  envelope.  This  spatial  pulse  was  seen  to 
require  about  txi70  orders  of  magnitude  improvement  in  the 
alignment  accuracy  of  the  optical  axes  of  the  lasers,  i.e.  on 
the  order  of  10~®  rad.  By  the  use  of  z transforms,  the  side 
lobes  of  the  spatial  pulse  could  be  controlled  in  any  desired 
manner.  For  the  monochromatic  case  (Af=0) , the  spatial  beam 
could  be  made  to  scan  sinusoidally  and  a planar  array  was 
seen  to  permit  conical  scanning.  These  results  make  the 
system  applicable  to  target  tracking  and  even  track-while- 
scan  since  the  beams'  relative  phases  can  bo  changed 
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electronically.  For  some  applications,  however,  the  second- 
ary spatial  beams,  called  grating  lobes,  are  undesirable,  so 
' the  cases  of  aperiodic  and  random  arrays  were  considered 

briefly.  Of  the  aperiodic  arrays,  the  spatial  taper  was 

found  to  provide  some  of  the  lowest  side  lobes,  but  for  op- 

j: 

tical  wavelengths  the  beam  pa'Jtern  away  from  the  main  peak 
became  less  controllable.  The  class  of  arrays  with  element 
locations  determined  in  a random  manner  from  array  to  array 
was  found  to  have  a ratio  of  average  side  lobe  intensity  to 
peak  intensity  inversely  proportional  to  the  number  of 
lasers.  Furthermore,  to  obtain  a reasonably  high  probability 
of  low  side  lobes  (i.e.  no  grating  lobes),  the  array  was 
found  to  require  a large  number  of  lasers  (>1000).  Although 
this  was  noted  to  be  true  in  general  for  the  class  of  random 
arrays,  the  results  only  guaranteed  that  a desirable  array 
existed  and  did  not  in  any  way  specify  the  array. 

Recommendations 

It  is  very  difficult  at  this  point  to  make  specific 
recommendations  as  to  what  type  of  system  should  be  built, 

, I how  many  lasers  it  should  use,  what  frequencies  it  should 

' I 

■ operate  at,  and  so  on.  The  difficulty  arises  because  as  is 

I 

{ usually  the  case,  none  of  the  issues  are  clear  enough  to 

) 

enable  instant  conclusions  to  be  drawn.  As  always,  there 
are  trade-offs  necessary  between  system  requirements,  design 
constraints,  and  costs.  The  difficulty  becomes  even  mere 
' acute  when  the  applications  are  generalized  as  they  are  here, 

i 

' I 
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and  when  there  are  no  systems  against  which  point  by  point 
comparisons  may  be  made.  Nevertheless,  some  general  comments 
may  be  made.  The  statistical  models  used  indicate  that 
given  a good  phase  locked  loop  design  for  heterodyne  detec- 
tion, the  fields  may  be  combined  and  detected  coherently. 
Misalignment  problems  in  the  plane  of  the  array  and  in  optic 
axis  angle,  though  difficult,  also  seem  to  be  within  the 
realm  of  engineering  solutions.  The  waist  alignment  problem 
(in  the  direction  of  propagation)  is  much  more  serious  and 
may  represent  the  most  serioua  stumbling  block  to  implement- 
ing coherent  combination  of  two  or  more  lasers.  The  flexi- 
bility of  such  a system  is  desirable  for  ranging  applications 
and  possibly  for  velocity  measurements  as  well  since  the 
pulse  format  may  be  changed  easily.  This  would  indicate, 
however,  that  more  than  two  or  a few  lasers  be  used  (exactly 
how  many  more  depends  on  the  system  requirements)  since  the 
measurement  precision  for  coupled  measurements  may  be  more 
easily  manipulated  for  N>2.  The  coaxial  case  requires  some 
sort  of  combination  scheme,  however,  which  may  be  lossy 
as  well  as  expensive.  A linear  array  may  find  application  in 
cases  where  beam  scanning  or  pointing  is  necessary.  However, 
unless  it  is  possible  to  build  arrays  of  more  than  several 
hundred  lasers  within  a length  of  a few  meters  or  less,  the 
multiple  lobe  problem  cannot  be  eliminated.  This  will  cer- 
tainly restrict  the  application  of  such  a system  even  though 
the  high  resolution  of  the  coherent  result  ( ~ cm)  may  be 
desirable . 
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There  are  several  areas  in  which  further  investigation 
would  be  helpful  and  in  which  hardware  development  should  be 
emphasized.  The  duality  issue  should  be  considered  more 
carefully.  It  appears  that  not  only  could  more  results  be 
obtained,  but  a great  deal  of  spatial  statistics  from  antenna 
theory  could  be  applied  to  the  temporal  realm.  As  far  as  the 
spatial  beam  is  concerned,  there  should  be  further  study  of 
the  spatial  taper  case  for  the  optical  wavelengths.  Either 
a general  theory  could  be  developed  or  the  approach  could  be 
shown  conclusively  to  be  of  little  value.  Computers  would 
be  helpful  in  determining  actual  patterns  for  various  arrays 
to  see  if  there  are  any  useful  trends.  For  the  random  case, 
the  statistics  should  be  examined  more  closely  to  produce 
more  rigorous  results.  Another  interesting  possibility  is 
that  of  developing  a measure  of  the  probability  that  an 
objectionably  high  side  lobe  will  occur  within  a certain 
region  near  the  main  beam.  Another  area  is  that  of  signal 
design.  Some  of  the  performance  criteria  presented  in  this 
study  are  basic  and  could  be  expanded.  Signal  synthesis 
would  also  open  up  the  whole  area  of  applications  to  communi- 
cations that  were  not  considered  here.  It  seems  .that  the 
coherent  combination  of  lasers  would  have  enough  more  flexi- 
bility than  a mode-locked  laser  to  be  useful  in  communica- 
tions. The  spatial  scanning  aspects  of  the  array  could  also 
be  helpful  in,  perhaps,  selecting  various  receivers  for 
various  messages,  i.e.  in  addition  to  frequency  and  time 
division  multiplexing,  spatial  m\jltiplexing  could  be  useful. 


Coherent  combination  of  lasers,  although  a formidable 
engineering  task,  still  offers  some  attractive  benefits. 
Development  of  arrays  of  small  lasers  with  power  outputs  of 
watts  or  milliwatts  should  be  encouraged.  Laser  alignment 
is  enough  of  a problem  that  it  is  desirable  to  obtain  more 
lasers  and  more  power  per  laser  in  a given  array.  Other 
hardware  problems  that  must  be  addressed  include  the  physical 
aspects  of  putting  a large  number  of  lasers  together  in  a 
small  array  with  the  associated  cooling,  control  electronics, 
and  optics. 
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